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Abstract 

^ ■ We derive the equations of motion for a planar rigid body of circular shape 

' moving in a 2D perfect fluid with point vortices using symplectic reduction by 

stages. After formulating the theory as a mechanical system on a configuration 
space which is the product of a space of embeddings and the special Euclidian 
group in two dimensions, we divide out by the particle relabelling symmetry and 
then by the residual rotational and translational symmetry. The result of the 
first stage reduction is that the system is described by a non-standard magnetic 
I symplectic form encoding the effects of the fluid, while at the second stage, a 

careful analysis of the momentum map shows the existence of two equivalent 
Poisson structures for this problem. For the solid-fluid system, we hence recover 
. the ad hoc Poisson structures calculated by Shashikanth, Marsden, Burdick and 

5— ( I Kelly on the one hand, and Borisov, Mamaev, and Ramodanov on the other 

^ ■ hand. As a side result, we obtain a convenient expression for the symplectic 

leaves of the reduced system and we shed further light on the interplay between 
curvatures and cocycles in the description of the dynamics. 
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1 Introduction 

In this paper, we use symplectic reduction to derive the equations of motion for a 
rigid body moving in a two-dimensional fluid with point vortices. Despite the fact 
that this setup is easily described, it may come as a surprise that these equations 

were derived (using straightforward calculations) only recently by Shashikanth, Marsden, Burdick, and Kell 
[2002] and Borisov, Mamaev, and Ramodanov [2003]. In order to see why this is so, 
despite the long history of this kind of problems, and in order to set the stage for our 
approach, let us trace some of the history of fluid-rigid body interaction problems. 

Rigid Bodies in Potential Flow. The equations of motion for a rigid body of 
mass TTib and inertia tensor 1;, were first described by Kirchhoff [1877] and are given 
by 

L = A X V and A = A x fi, (1.1) 



1 Introduction 3 



where V and Q are the Unear and angular velocity of the body, while L and A 
are the linear and angular momentum, related by L = mV and A = Ift. Here, 
m = rrib + m/ is the total mass of the rigid body, consisting of the mass rrih and the 
virtual mass mj induced by the fluid. Similarly, I = 1;, + 1^, where 1/ is the virtual 
inertia tensor due to the fluid. 

The main difference between these equations and the Euler equations governing the 
motion of a rigid body in vacuum is the appearance of the non-zero term A x V 
on the right-hand side of the equation for L. In other words, the center of mass no 
longer describes a uniform straight trajectory and is a non-trivial degree of freedom. 
From a geometric point of view, the motion of a rigid body in a potential flow can 
therefore be considered as a curve on the special Euclidian group SE(3) consisting 
of all translations and rotations in M^, where the former describe the orientation of 
the body while the latter encode the location of the center of mass. 

The kinetic energy for the rigid body in a potential flow is a quadratic function 
and hence determines a metric on SE(3). The physical motions of the rigid body 
are geodesies with respect to this metric. By noting that the dynamics is invariant 
under the action of SE(3) on itself, the system can be reduced from the phase space 
T* SE(3) down to one on the dual Lie algebra 5e(3)*: in this way, one derives the 
Kirchhoff equations (1.1). The situation is similar for a planar rigid body moving in 
a two-dimensional flow: the motion is then a geodesic flow on SE(2), the Euclidian 
group of the plane, and reduces to a dynamical system on 5e(2)*. The dynamics 
of the reduced system is still described by the Kirchhoff equations, but additional 
simplifications occur: since both A and Ct are directed along the z-axis, the right- 
hand side of the equation for A is zero. Throughout this paper, we will deal with 
the case of planar bodies and two-dimensional flows only. 

The procedure of reducing a mechanical system on a Lie group G whose dynamics 
are invariant under the action of G on itself is known as Lie-Poisson reduction. 
It was pointed out by Arnold [1966] that both the Euler equations for a rigid body 
as well as the Euler equations for a perfect fluid can be derived using this approach. 
For more about the history of these and related reduction procedures, we refer 
to Marsden and Ratiu [1994]. The geometric outlook on the Kirchhoff equations, 
developed by Leonard [1997] , turned out to be crucial in the study stability results 
for bottom-heavy underwater vehicles; see also Patrick, Roberts, and Wulff [2008]. 

Point Vortices. Another development in fluid dynamics to which the name of 
Kirchhoff is associated, is the motion of point vortices in an inviscid flow. A point 
vortex is defined as a singularity in the vorticity field of a two-dimensional fiow: 
Lo = r(5(x — xo(t)), where the constant F is referred to as the strength of the point 
vortex. By plugging a superposition of N point vortices into the Euler equations 
for a perfect fluid, one obtains the following set of ODEs for the evolution of the 
vortex locations Xj, i = 1, . . . , iV: 



= J 



dH 



(1.2) 




9xfe' 
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where H = — VFg'(xi, . . . , xjv) and Wq is the so-cahed Kirchhoff-Routh func- 
tion, which is derived using the Green's function for the Laplacian with appropriate 
boundary conditions. We wih return to the exphcit form for Wg later. For point 
vortices moving in an unbounded domain, Wg was derived by Routh and Kirchhoff, 
whereas Lin [1941] studied the case of vortices in a bounded container with fixed 
boundaries. 

It is useful to recall here how the point vortex system is related to the dynamics of 
an inviscid fluid. In order to do so, we recall the deep insight of Arnold [1966], who 
recognized that the motion of a perfect fluid in a container ^ is a geodesic on the 
group Diffvoi(-7^) of volume-preserving diffeomorphisms of J-', in much the same way 
as the motion of a rigid body is a geodesic on the rotation group S0(3). The group 
Diffvoi(.^) acts on itself from the right and leaves the fluid kinetic energy invariant, a 
result known as particle relabelling symmetry. Hence, by Lie-Poisson reduction, 
the system can be reduced to one on the dual of the Lie algebra of Diffvoi(-^), 

and the resulting equations are precisely Euler's equations. 

Moreover, any group acts on its dual Lie algebra through the co-adjoint action, and a 
trajectory of the reduced mechanical system starting on one particular orbit of that 
action is constrained to remain on that particular orbit. Marsden and Weinstein 
[1983] showed that for a perfect fluid, the co-adjoint orbits are labelled by vorticity 
and, when specified to the co-adjoint orbit corresponding to the vorticity of N point 
vortices, Euler's equations become precisely the vortex equations (1.2). 

For a detailed overview of the geometric approach to fluid dynamics, we refer to 
Arnold and Khesin [1998]. The dynamics of point vortices is treated by different 
authors, most notable Saffman [1992], Newton [2001] and Aref [2007]. 

The Rigid Body Interacting with Point Vortices. Given the interest in 
both rigid bodies and point vortices, it is natural to study the dynamics of a 
planar rigid body interacting with N point vortices. Surprisingly, it wasn't until 
the recent work of Shashikanth, Marsden, Burdick, and Kelly [2002] (SMBK) and 
Borisov, Mamaev, and Ramodanov [2003] (BMR) that the equations of motion for 
this dynamical system were established. Both groups proceeded through an ad-hoc 
calculation to derive the equations of motion, and showed that the resulting equa- 
tions are in fact Hamiltonian. However, both sets of equations are different at first 
sight. 

The SMBK equations are formally identical to the Kirchhoff equations (1.1) together 
with the point vortex equations (1.2), but the definitions of the momenta L, A and 
the Hamiltonian H are modified to include the effect of the rigid body (through 
the ambient fluid) on the point vortices, and vice versa. The configuration space 
is se(2)* X and the Poisson structure is the sum of the Poisson structures on 
the individual factors. From the BMR point of view, the dynamic variables are the 
velocity V and angular velocity A, together with the locations of the vortices. The 
Hamiltonian is simply the sum of the kinetic energies for both subsystems, and to 
account for the interaction between the point vortices and the rigid body, BMR 
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introduce a non-standard Poisson structure involving the stream functions of the 
fluid. 

Somewhat miraculously, the equations of motion obtained by SMBK and BMR 
turn out to be equivalent: Shashikanth [2005] establishes the existence of a Poisson 
map taking the canonical SMBK Poisson structure into the BMR Poisson structure. 
However, a number of questions therefore remain. Most importantly, it is not ob- 
vious why one mechanical problem should be governed by two Poisson structures, 
which are at first sight very different but turn out to be related by a certain Poisson 
map. Moreover, it is entirely non-obvious why this system is Hamiltonian in the 
first place: in the work of SMBK and BMR, the Hamiltonian structure is derived 
only afterwards by direct inspection. 

Main Contributions and Outline of this Paper. We shed more light on the 
issues addressed above by uncovering the geometric structures that govern this prob- 
lem. We begin this paper by giving an overview of rigid body dynamics and aspects 
of fluid mechanics in section 2. This material is mostly well-known and serves to 
set the tone for the rest of the paper. In section 3 we then introduce the so-called 
Neumann connection, giving the response of the fluid to an infinitesimal motion of 
the rigid body. This connection has been described before, but a detailed overview 
of its properties seems to be lacking. In particular, we derive an expression for 
the curvature of this connection when the fiuid space is an arbitrary Riemannian 
manifold, generalizing a result of Montgomery [1986]. 

The remainder of the paper is then devoted to using reduction theory to obtain the 
equations of motion for the fluid-solid system in a systematic way. We will derive 
the equations of motion for this system by reformulating it first as a geodesic flow on 
the Cartesian product of the group SE(2) of translations and rotations in 2D, and 
a space of embeddings describing the fluid configurations. Two symmetry groups 
act on this space: the particle relabelling symmetry group Diffvoi(-^), and the group 
SE(2) itself. Dividing out by the combined action of these symmetries is therefore 
an example of reduction by stages. 

After symplectic reduction with respect to DiffyoK^) in section 4, we obtain a 
dynamical system on the product space T* SE(2) x M."^^ , where the first factor is 
the phase space of the rigid body, while the second factor describes the locations of 
the point vortices. The dynamics is governed by a magnetic symplectic form: 
it is the sum of the canonical symplectic forms on both factors, together with a 
non-canonical interaction term. From a physical point of view, the interaction term 
embodies the interaction between the point vortices and the rigid body through the 
ambient fluid. Mathematically speaking, the map associating to each rigid body 
motion the corresponding motion of the fluid can be viewed as a connection, and 
the magnetic term is closely related to its curvature. 

Note that it is essential here to do symplectic reduction with respect to the particle 
relabeling symmetry rather than Poisson reduction, even though the latter might 
be conceptually simpler. Recall that symplectic reduction reduces the system to a 
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Figure 1.1: The two symmetry groups in the problem of a cyhnder interacting with A'' point 
vortices are the particle relabelling group DiS^oii^) and the Euclidian group SE(2). 

CO- adjoint orbit; as we shall see below, in the case of fluids interacting with solids 
these co-adjoint orbits are precisely labelled by the vorticity of the fluid. Symplectic 
reduction — in particular the selection of one particular level set of the momentum 
map — therefore amounts to imposing a specific form for the vorticity field. In our 
case, this is where we introduce the assumption that the vorticity is concentrated 
in point vortices. In the case of Poisson reduction, we would have obtained the 
equations of motion for a rigid body interacting with an arbitrary vorticity field. 

After factoring out the particle relabeling symmetry, the resulting dynamical system 
is invariant under translations and rotations in the plane and can then be reduced 
with respect to the group SE(2). This is the subject of section 5. Physically, this 
corresponds to rewriting the equations of motion obtained after the first reduction 
in body coordinates. However, because of the presence of the magnetic term in the 
symplectic form and the fact that SE(2) acts diagonally on T* SE(2) x M^^, this is 
not a straightforward task. 

First, we derive the reduced Poisson structure on the reduced space 5e(2)* x M?^ . 
Because of the magnetic contributions to the symplectic form, the reduced Poisson 
structure is not just the sum of the Poisson structures on the individual factors, 
but includes certain non-canonical contributions as well. We then show that the 
momentum map for the SE(2)-symmetry naturally defines a Poisson map taking 
this Poisson structure into the product Poisson structure, possibly with a cocycle if 
the momentum map is not equivariant (this happens when the total strength of the 
point vortices is nonzero). We do the computations for a general product T*H x P 
first, where H isa Lie group, P is a symplectic manifold, and the product is equipped 
with a magnetic symplectic form. In this way, we generalize the "coupling to a Lie 
group" scenario (see Krishnaprasad and Marsden [1987]) to the case where magnetic 
terms are present. 

In this way, the results of SMBK and BMR are put on a firm geometric footing: 
the BMR Poisson structure is the one obtained through reduction and involves the 
interaction terms, while the SMBK Poisson structure is simply the product Poisson 
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structure. The Poisson map induced by the SE(2)-momentum map described above 
turns out to be precisely Shashikanth's Poisson map. As a consequence, we also 
obtain an explicit prescription for the symplectic leaves of this system. 

Relation with Other Approaches. Our method consists of rederiving the SMBK 
and BMR equations by reformulating the motion of a rigid body in a fluid as a 
geodesic problem on the space Q. By imposing the assumption that the vorticity 
is concentrated in N point vortices, and dividing out the symmetry, we obtain first 
of all the BMR equations, and secondly (after doing a momentum shift) the SMBK 
equations. This procedure is worked out in the body of the paper — here we would 
like to reflect on the similarities with other dynamical systems. 

Recall that the dynamics of a particle of charge e in a magnetic field B = V x A 
can be described in two ways. One is by using canonical variables (q, p) and the 
Hamiltonian 

^(q,p) = ;^l|p-eAf , (1.3) 

while for the other one we use the kinetic energy Hamiltonian -ffkin = ||p||^ /2m but 
now we modify the symplectic form to be 

f^B = ^^can + eB. (1.4) 

A simple calculation shows that both systems ultimately give rise to the familiar 
Lorentz equations. From a mathematical point of view, this can be seen by noting 
that the momentum shift map 5 : p i— > p — eA maps the dynamics of the former 
system into that of the latter: 

5* Ob = Ocan and 5* //kin = H. 

In other words, both formulations are related by a symplectic isomorphism, thus 
making them equivalent. 

An over-arching way of looking at the dynamics of a charged particle in a mag- 
netic field is as a geodesic problem (with respect to a certain metric) on a higher- 
dimensional space: this is part of the famed Kaluza-Klein approach. In this 
case, spacetime is replaced by the product manifold x U(l) and charged parti- 
cles trace out geodesies on this augmented space. The standard, four-dimensional 
formulation of the dynamics can then be obtained by dividing out by the internal 
U(l)-symmetry, resulting in the familiar Lorentz equations on M^. More information 
about these constructions can be found in Sternberg [1977], Weinstein [1978] and 
Marsden and Ratiu [1994] . 

Our approach to the fluid-structure problem is similar but more involved because of 
the presence of two different symmetry groups. However, the underlying philosophy 
is the same: by reformulating the motion of a rigid body and the fluid as a geodesic 
problem on an infinite-dimensional manifold, we follow the philosophy of Kaluza- 
Klein of trading in the complexities in the equations of motion for an increase in the 
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number of dimensions of the configuration space. Then, by dividing out the Diffyoi- 
symmetry, we obtain a reduced dynamical system governed by a magnetic symplectic 
form (the analogue of (1-4)), which is mapped after a suitable momentum shift S 
into an equivalent dynamical system governed by the canonical symplectic form but 
with a modified Hamiltonian. After reducing by the residual SE(2)-symmetry, the 
former gives rise to the BMR bracket, while the latter is nothing but the SMBK 
system. 

A Note on Integrability. The case of a circular cylinder is distinguished because 
of the existence of an additional conservation law, associated with the symmetry 
that rotates the cylinder around its axis. When the circular cylinder interacts with 
one external point vortex, a simple count of dimensions and first integrals suggests 
that this problem is integrable, a fact first proved by Borisov and Mamaev [2003]. 
Indeed, for a single vortex the phase space is a symplectic leaf of se*(2) x M^, which 
is generically four-dimensional. On the other hand, for the circular cylinder there 
are two conservation laws: the total energy and the material symmetry associated 
with rotations around the axis of the cylinder, which hints at Liouville integrability. 

In the case of an ellipsoidal cylinder with nonzero eccentricity, Borisov, Mamaev, and Ramodanov 
[2007] gather numerical evidence to show that the interaction with one vortex is 
chaotic. This is to be contrasted with the motion of point vortices in an unbounded 
domain (see Newton [2001]), which is integrable for three vortices or less. 
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2 The Fluid-Solid Problem 

This section is subdivided into two parts. In the first, and longest, part we describe 
the general setup for a planar rigid body interacting with a 2D flow. The second 
part is then devoted to discussing a number of simplifying assumptions that will 
make the subsequent developments clearer. By separating these assumptions from 
the main problem setting, we hope to convince the reader that the method outlined 
in this paper does not depend on any specific assumptions on the rigid body or the 
fluid, and can be generalized to more complex problems. At the end of this paper. 
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we discuss how these assumptions can be relaxed. 
2.1 General Geometric Setting 

Kinematics of a Rigid Body. Throughout this paper, we consider the motion 
of a planar rigid body interacting with a 2D flow. We introduce an inertial 
frame ei.2,3, where ei_2 span the plane of motion and 63 is perpendicular to it. The 
configuration of the rigid body is then described by a rotation with angle /3 around 
e3 and a vector xq = x^ei +2/0^2 describing the location of a fixed point of the body, 
which we take to be the center of mass. 




X 



Figure 2.1: The configuration of the rigid body and the point vortices. 

The orientation and position (/3;a;o,yo) determine an element g of the Euclidian 
group SE(2) given by 

5=fn where i?=f^°^^ " V (2.1) 

\0 I J \sm(3 cos (3 J ^ ^ 

Written in this way, the group composition and inversion in SE(2) are given by 
matrix multiplication and inversion. The Lie algebra of SE(2) is denoted by 5e(2) 
and essentially consists of infinitesimal rotations and translations. Its elements are 
matrices of the form 

-n vA 

n Vy\ . 

00 / 

It follows that se(2) can be identified with M'^ by mapping such a matrix to the 
triple (O, Vx, Vy). We define a basis of se(2) given by 

/O -1 0\ /O 1\ /O 0\ 

ef^= l , e^= , 6^= 1. (2.2) 
\0 0/ VO 0/ \0 0/ 

For future reference, we also introduce a moving frame fixed to the rigid body, 
denoted by bi^2,3- The transformation from body to inertial frame is given by 
X = i?X + xq, where x = xei + ye2 and X = Xhi + Yh2- 
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The angular and translational velocities of the rigid body relative to the inertial 
frame are defined as 

uj = /3e3 and v = Xoei + 1)062, (2.3) 

where dots denote derivatives with respect to time. These quantities can be ex- 
pressed in the body frame: the body angular velocity fl and the body translational 
velocity V are related to the corresponding inertial quantities by 

ft = u and V = R^v. (2.4) 

For the case of a planar rigid body, the angular velocity is oriented along the axis 
perpendicular to the plane and so determines a scalar quantity: lj = u>e^ and 
CI = Oba. From a group-theoretic point of view, if the motion of a rigid body is 
given by a curve t 1— > g{t) in SE(2), then we may define an element (,{t) € se{2) by 
putting 

U,€ e=U 

and it can easily be checked that ^ coincides with the body angular and translational 
velocities (r2,V). 

The kinetic energy of the rigid body is given by the following expression: 

rbody(^^,V) = ^J72 + ^V2 (2.5) 

where I is the moment of inertia of the body and m is its mass. The kinetic energy 
defines an inner product on se(2), given by 

(((Ol, Vi), (^^2, V2))),,(2) = Vi) Mm (^'J , (2.6) 

for (Oi,Vi) and (J^2,V2) in se(2), and where M^, is given by 

1 



mlj ■ 

Here, I is the 2-by-2 identity matrix. By left extension, this inner product induces 
a left-invariant metric on the whole of SE(2): 

(((5',5'i),(fi',5'2)))sE(2) ■= ((5'~^5i,5"^5'2»g,(2) , (2-8) 
for {g,gi) and {9,92) in rSE(2). 



Incompressible Fluid Dynamics. The geometric description of an incompress- 
ible fluid goes back at least to the work of Arnold [1966], who described the motion 
of an incompressible fluid in a fixed container ^ as a geodesic on the diffeomorphism 
group of J^. In the case of a fluid interacting with a rigid body, the fluid container 
may change over time, reflecting the fact that the rigid body moves. 
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Arnold's formulation can be extended to cover this case by considering the space 
of embeddings of the reference configuration of the fluid (denoted by ^o) into M^. 
Recall that an embedding ip : J'q ^ maps each reference point X G J^q to its 
current configuration x = (p{X). In order to reflect the fact that the fluid is taken to 
be incompressible, we require that any fluid embedding ip is volume preserving: 
if dx is the Euclidian area element on M? and i]q is a fixed area element on !Fq, then 
we require 

V9*(dx) = r?o. (2.9) 

The space of all such volume-preserving embeddings is denoted by Emb(.7-o, M^). In 
the sequel, we will specify additional boundary conditions on the fluid configurations, 
stating for example that the fluid is free to slide along the boundary of the solid, 
but this does not make any difference for the current expository treatment. 

A motion of a fluid is described by a curve t i-^ ip{t) in Fjmh{!FQ,M?). The material 
velocity field is the tangent vector field {<pt, ift) along the curve. Here, ipt is a map 
from To to TR^, whose value at a point x £ J^o is given by 

Note that (pt is not a vector field in the traditional sense; rather, it is a vector field 
along the map ipt- In contrast, the spatial velocity field u^, defined as 

Ut = (ptOip~'^ (2.10) 

is a proper vector field, defined on Lpt{J^Q). 

The motion of a fiuid can be described using the kinetic-energy Lagrangian: 

Tf,uid{'P,'P) = ^ l|u|p(ix, (2.11) 

where u is the Eulerian velocity field (2.10) and the integration domain is the spatial 
domain of the fiuid at time t: Tt = '■Pti^o)- Just as in the case of the rigid body, 
this kinetic energy induces a metric on the space Emb(.7^0) I^^)) given by 

((VJl, '/j2))Emb = P 'Pl-'P2 m- (2-12) 

By changing variables, this metric can be rewritten in spatial form as follows: 

(('^l,'/i2))Emb = P I • Uadx, 

where Uj is the Eulerian velocity associated to (/9j: Uj = ipi o ip^^ , for i = 1,2, and 
the integration is again over the spatial domain J^t ■= Vti^o) of the fluid. 
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The Configuration Space of the Fluid-Solid System. The motion of a rigid 
body in an incompressible fluid combines aspects of both rigid-body and fluid dy- 
namics. We assume that the body occupies a circular region Bq in the reference 
configuration, and that the remainder of the domain, denoted by JFq, is taken up 
by the fluid. The configuration space for the fluid-solid system is made up of pairs 
{g,ip) G SE(2) X Emb(jFo,]R^) satisfying the following conditions. 



1. The embedding (p represents the configuration of the fluid. In particular, (p 
is volume-preserving, i.e. (2.9) is satisfled. In addition, we assume that ip 
approaches the identity at inflnity suitably fast.^ 

2. The element g G SE(2) describes the conflguration of the rigid body. 

3. The fluid satisfles a "slip" boundary condition: the normal velocity of the 
fluid coincides with the normal velocity of the solid, while the tangent ve- 
locity can be arbitrary, reflecting the fact that there is no viscosity in the 
fluid. Mathematically, this boundary condition is imposed by requiring that, 
as sets, ip{dTo) is equal to g{dI3o), where g € SE(2) is interpreted as a linear 
embedding of Bq into M^. 

We denote the space of all such pairs as Q; this is a submanifold of SE(2) x 
Emb(J^o,]R^)- The kinetic energy of the fluid-solid system is given by the sum 
of the rigid-body energy Tbody and the kinetic energy Tfjuid of the fluid: 

T = Tfluid + Tbody 

= ^ [ \\nfd^ + -n' + -V^ (2.13) 
2^" " 2 2 ^ ' 

Similarly, there exists a metric on Q given by the sum of the metrics (2.8) and 
(2.12): 

{{{^l,9l),{^2,g2))) = ((V^l > 9^2 ))Emb + ((51,52)) SE(2) ' (2-14) 

The dynamics of rigid bodies moving in perfect fluids was studied before by Kelly 
[1998], Radford [2003], and Kanso, Marsden, Rowley, and Melli-Huber [2005]; Kanso and Oskouei 
[2008]. A similar conflguration space, but with the SE(2)-factor replaced by a suit- 
able set of smooth manifolds, was studied by Lewis, Marsden, Montgomery, and Ratiu 
[1986] for the dynamics of a liquid drop. 



Particle Relabelling Symmetry. The kinetic energy Tfjuid of the fluid is invari- 
ant if we replace p hy (p o (p, where (f> is a volume-preserving diffeomorphism from 
J-'o to itself. This represents the particle relabelling symmetry referred to in the 
introduction. 

^We will not be concerned with any functional-analytic issues concerning these infinite- 
dimensional manifolds of mappings. Instead, the reader is referred to Ebin and Marsden [1970] 
for more information. 
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Recall that a diffeomorphism (p : ^ Tq \s volume-preserving if = where 
r/o is the volume element on Tq. The group of all volume-preserving diffeomorphisms 
is denoted by Diffvoi- This group acts on the right on Emb(^05l^^) by putting 
ip ■ 4> = (f o and hence also on Q. The action of Diffvoi on Q makes Q into the 
total space of a principal fiber bundle over SE(2). In other words, if we define the 
projection tt : Q —i- SE(2) as being the projection onto the first factor: Tr{g, ip) = g, 
then the fibers of vr coincide precisely with the orbits of Diffvoi in Q- 



Vorticity and Circulation. In classical fluid dynamics, the vorticity n is de- 
fined as the curl of the velocity field: = V x u, and the circulation around the 
rigid body is the line integral of u along any curve encircling the rigid body. In two 
dimensions, fi can be written as = fie^, where /i is called the scalar vorticity. 

According to Noether's theorem, there is a conserved quantity associated to the 
particle relabelling symmetry. This conserved quantity turns out to be precisely the 
circulation of the fluid: 

— i u • dl = 0, 

dtJc 

and Noether's theorem hence becomes Kelvin's theorem, which states that circula- 
tion is materially constant. As a consequence of Green's theorem, the circulation of 
the fluid is related to the vorticity: 



u dl= / V X u • dS, 

c Js 

where S is any surface whose boundary is C. It would lead us too far to explore the 
geometry of vorticity and circulation in detail (for this, we refer to Arnold and Khesin 
[1998]) but at this stage we just note that the conservation of vorticity is closely 
linked to the particle relabelling symmetry. 



The Helmholtz-Hodge Decomposition. If (in addition to being incompress- 
ible) the fluid is irrotational, meaning that V x u = 0, then there exists a velocity 
potential $ such that u = V$. In the presence of point vortices, the fluid is not 
irrotational but the velocity field u can be uniquely decomposed in an irrotational 
part and a vector field Uv representing the "rotational" contributions (see for in- 
stance Saffman [1992] or Newton [2001]). This is the well-known Helmholtz-Hodge 
decomposition: 

u = V«> + Uv. (2.15) 

Here, Uv is a divergence-free vector field which is tangent to the boundary of J^, 
while the potential ^ is the solution to Laplace's equation V^<I> = 0, subject to the 
boundary conditions that the normal derivative of $ equals the normal velocity of 
the rigid body, and that the velocity vanishes at infinity. In other words, one has 



= and — 
on 



= (cj X (x-xo) + v) -n (foTxedT), (2.16) 
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and $ goes to zero as |x| goes to infinity. Here u and v are the angular and trans- 
lational velocity of the rigid body, expressed in a spatial frame, while xq represents 
the location of the center of mass. 

Since <1> depends on lo and v in a linear way, we may decompose following 
Kirchhoff, as 

^>(5,g;x) = v^$^(c/,x) +7;j/$j^(c/,x) +w^>^(c/,x). (2.17) 

Here ^y, and $ are elementary velocity potentials corresponding to infinites- 
imal translations in the x- and y-direction and to a rotation, respectively. They 
satisfy the Laplace equation with the following boundary conditions: 

— — = (e3 X x) • n, — — = ei • n, and — = 62 • n. 

on on on 

Note that ^ depends on the boundary data, and hence on (g, g), while the elementary 
potentials depend on the location of the rigid body (encoded by g) only. For the 
sake of clarity, we will suppress these arguments when no confusion is possible. The 
elementary velocity potentials for a circular body are given below: the reader can 
check that they indeed depend on the location of the rigid body. 

The Helmholtz-Hodge decomposition defines a connection on the principal fiber 
bundle Q — > SE(2). To see this, note that any tangent vector to Q is of the form 
(g,g;(p,ip). We define its horizontal and vertical part by applying the Helmholtz- 
Hodge decomposition to u = c/i o ip~^, and we put 

{g,g]V,^)" = {9,9;^^^^ and {g,g]ip,ifY = {g,Q;Lp,u^ o ^p) 

We will verify in section 3 below that this prescription indeed defines a connection, 

which we term the Neumann connection, since its horizontal subspaces are found 

by solving the Neumann problem (2.16). This connection was used in a variety of 

contexts, ranging from the dynamics of fluid drops (see Lewis, Marsden, Montgomery, and Ratiu 

[1986]; Montgomery [2002]) to problems in optimal transport (see Khesin and Lee 

[2008]). 



The Lie Algebra of Divergence-free Vector Fields and its Dual. At least 
on a formal level, DifFvoi is a Lie group with associated to it a Lie algebra, denoted by 
3£voi and consisting of divergence-free vector fields which are parallel to the boundary 
of J-Q. The bracket on Xvoi is the Jacob i-Lie bracket of vector fields (which is the 
negative of the usual bracket of vector fields) and its dual, denoted by X*^!, is the 
set of linear functionals on Xvoi- This set can be identified with the set of one- forms 
on T{) modulo exact forms: 

x:„i(J-o) = Oi(.Fo)/df]0(.Fo), 

(see Arnold and Khesin [1998]) while the duality pairing between elements of Xvoi 
and is given by 

(u, [a]) = / a(u)r/o. 
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where u G X^oi and [a] G ^voi- Note that the right-hand side does not depend on 
the choice of representative a. 

Another interpretation of the dual Lie algebra X*^i is as the set of exact two-forms 
on Tq. Any class [ Q^] in -^^Q^ is uniquely determined by the exterior differential dct 
and by the value of a on the generators of the first homology of J-q. In our case, the 
first homology group of jFg is generated by any closed curve C encircling the rigid 
body, and its pairing with a is given by 



Under this identification, da represents the vorticity , while T represents the cir- 
culation. Since in our case the circulation is assumed to be zero, it follows that [a] 
is completely determined by da. 

From a geometric point of view, the vorticity field can be interpreted as an element 
of if we assume that the fluid is moving on an arbitrary Riemannian manifold, 
then the vorticity can be defined by 



Here b is the flat operator associated to the metric. In the case of Euclidian spaces, 
this definition reduces to the one involving the curl of u. This definition is slightly 
different from the one in Arnold and Khesin [1998], where vorticity is defined as 
a two- form on the inertial space M?, whereas in our interpretation, vorticity lives 
on the material space J-'q. Both definitions are related by push- forward and pull- 
back by (f, and hence carry the same amount of information. Our definition has 
the advantage that fi is naturally an element of X*^^, which is preferable from a 
geometric point of view. 

We finish this section by noting that any Lie group acts on its Lie algebra and its 
dual Lie algebra through the adjoint and the co-adjoint action, respectively. For the 
group of volume-preserving diffeomorphisms, both are given by pull-back: if (p is an 
element of Diffvoi, and u and da are elements of jCvoI and X*^^, respectively, then 



2.2 Point Vortices Interacting with a Circular Cylinder 

In this section, we impose some specific assumptions on the rigid body and the 
fluid. It should be pointed out that while these assumptions greatly simplify the 
exposition, the general reduction procedure can be carried out under far less strin- 
gent assumptions. Later on, we will discuss how some of these assumptions may be 
removed. 

The Rigid Body. In order to tie in this work with previous research efforts, we 
assume the rigid body to be circular with radius R and neutrally buoyant {i.e. its 





Ad0(u) = (j)*u and CoAd^(da) = d(0*a). 
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body weight m is balanced by the force of buoyancy). If the density of the fluid is 
set to p = 1, this implies that the body has mass m = ttE?. The moment of inertia 
of the body around the axis of symmetry is denoted by I. 

For a rigid planar body of circular shape, the elementary velocity potentials 

and occurring in (2.17) can be calculated analytically (see Lamb [1945]) and are 

given by 

^- = -^7 wrf ^' % = -Ri Irzf (2-18) 

(x - xqY + {y - yor - + [y - yor 

while = 0, reflecting the rotational symmetry of the body. Here (2:0,2/0) the 
coordinates of the center of the disc. 

In some cases, it will be more convenient to express ^ in body coordinates. In 
analogy with (2.17), we may write 

where (V^, Vy) and O are the translational and angular velocity in the body frame, 
respectively. For the circular cylinder, the elementary potentials in body frame are 
given by 

^^ = -^^2^72' ^^ = -^x2^ and $n = 0. (2.19) 

Note that ^x,^Y and $q do not depend on the location of the rigid body, in 
contrast to ^x,^y and <I>^. 



Point Vortices. As for the fluid, we make the fundamental assumption that the 
vorticity is concentrated in N point vortices of strengths Fj, i = 1,. . . ,N, and 
that there is no circulation. Considered as a two- form on J^q, the former means that 
the vorticity is given by 

N 

I^ = Y1 ri5(x - 5ti)dx A dy, (2.20) 

i=l 

where {x, y) are coordinates on Tq, and Xj is the reference location of the ith vortex, 
i = 1, . . . , iV. As pointed out above, fi is an element of ^^qJ' 



The Kirchhoff-Routh Function. Shashikanth, Marsden, Burdick, and Kelly [2002] 
showed that the kinetic energy for the vortex system is the negative of the Kirchhoff- 
Routh function Wq for a system of N point vortices moving in a domain with moving 
boundaries: 

^vortex (Xi,...,X,v) = -W^G(Xi,...,XAr). (2.21) 
The precise form of Wg is given by 

Wg(Xi,...,X;v) = rir,G(x„x,) + i J^Fi<7(Xi,xo, (2.22) 
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where G(Xo,Xi) is a Green's function for the Laplace operator of the form 

1 2 
G(Xo,Xi) = 5f(Xo,Xi) + ^log ||Xo - Xill . 



(2.23) 



The function g is harmonic in the fluid domain and is the stream function of u/. As 
pointed out in Shashikanth, Marsden, Burdick, and Kelly [2002], in the case of a cir- 
cular cylinder, g can be calculated explicitly using Milne-Thomson's circle theorem 
(Milne-Thomson [1968]), and is given by 



5(X,Y) = i-log||Xf-^lo, 



X 



i?2 



From a geometrical point of view, Marsden and Weinstein [1983] showed that the 
vortex energy can be obtained (up to some "self-energy" terms) by inverting the 
relation = V x u, where // is the scalar part of (2.20), and substituting the 
resulting velocity field u generated by N vortices into the expression for the kinetic 
energy of the fluid. Although the analysis of Marsden and Weinstein [1983] was 
for an unbounded fluid domain, their result can easily be extended to the case 
considered here, yielding again the negative of the Kirchhoff-Routh function as in 
(2.21). 



Dynamics of the Fluid-Solid System. As stated in the introduction, Shashikanth, Marsden, Burdick, 

[2002] (SMBK) were the first to derive the equations of motion for a rigid cylinder 

interacting with N point vortices. These equations generalize both the Kirchhoff 

equations for a rigid body in a potential flow and the equations for point vortices 

in a bounded flow. Rather remarkably, SMBK established by direct inspection that 

these equations are Hamiltonian with respect to the canonical Poisson structure 

on se(2)* x M?^ [i.e. the sum of the Poisson structures on both factors), and a 

Hamiltonian given below involving the kinetic energy and interaction terms. 

The SMBK equations are given by 

-}r = 0, — + VxL = 0, and Pfc— ^ = -J-— , 2.24 
dt dt dt aXfc 

where L and A = Ae^ are the translational and angular momenta of the system, 
defined by 

N ^ X 

L = cV + ^rfcXfc xe3 + J^Ffceg X — ^ (2.25) 



k=l k=l 

N 



iXfci 



A = m-^^Ti\\x,f, 
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and H is the Hamiltonian: 

H{L,Xk) = -W{L,Xk) + ^ ||Lf - ^( J^rfe(L X Xfc) • eg 

-^Y^^l llXfcf - J^Tfcr.Xfc . X,- (2.26) 

+ 



j>k 

Xfc ^-^ „ Xfc 



XI I Z ' / ^ '" II ■\7- 1 1 2 



Here, c is the total mass of the cyhnder, consisting of the intrinsic mass m and the 
added mass ttR^ (due to the presence of the fluid): c = m + -kB? = 271 E? . 



Theorem 2.1. The SMBK equations are Hamiltonian on the space 5e(2) 
equipped with the Poisson bracket 



p2N 



{-^) G}sc{2)* xR2]V — {-F|sc(2)* ) G'|sc(2)* }l.p + {-f'lIRaiV , G^^2N jvortex- (2.27) 

Here, F and G are functions on se(2)* x M?^ . 

In the theorem above, {•, - ji p is the Lie-Poisson bracket on 5e(2)*: 

/ -Py PA 

{/i,/2W = (V/i)^AV/2, where A = , (2.28) 

\-Px / 

for arbitrary functions /i,/2 on 5e(2)*. Similarly, {-j-lvortex is the vortex bracket, 
given by: 

r . I f dgi dg2 dg2 dgi \ 

i=l 

where gi , g2 are arbitrary functions on . 

The BMR Equations. A completely different perspective on the rigid body 
interacting with N point vortices is offered by Borisov, Mamaev, and Ramodanov 
[2003] (BMR). From their point of view, the equations of motion are again written 
in Hamiltonian form on 5c(2)* x R^^, but now with a noncanonical Poisson bracket: 

F = {F,H}bmr 

for all functions F on 5e(2)* x M^^. The Hamiltonian is the sum of the kinetic 
energies of the subsystems, without interaction terms: 



H=^{Y,Y)-WG(V,Xk), 
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whereas the Poisson bracket is determined by its value on the coordinate functions: 

FRf-R^ lRf-R^{Xf-Y^) 



{Vl, V^2}bMR = -2 - X] 72 ' {Vi,Xi}BMR=^ 



Rf ' L'i,^^«jmv.^ ^ ^4 

{14,^i}BMR = {V2,-'^j}BMR = (2-30) 

ii/2, J^ijBMR — , {Xi,yiiBMR — -JT^ 

where i?^ = ||Xj|p, V = {Vi, V2), Xj = (Xj, yj), and T is the total vortex strength: 

AT 



I ■ 

i=l 



Note that this Poisson bracket differs from the one in Borisov, Mamaev, and Ramodanov 
[2003] by an overall factor of 27r. Shashikanth [2005] showed that this discrepancy 
can be attributed to the way in which BMR choose the fluid density. 

The Link between the SMBK and the BMR Equations. By explicit calcu- 
lation, Shashikanth [2005] constructed a Poisson map taking the SMBK equations 
into the BMR equations. His result is listed below. 

Theorem 2.2. The map S : (L, A; Xi, . . . , X^r) 1— > (V, fi; Xi, . . . , X/yr), where L 
and A are given by (2.25), is a Poisson map from se(2)* x equipped with the 
bracket {•, ■}0e(2)* xR^^v to se(2)* x with the bracket {•, -Ibmr- 

Even though this result asserts that both sets of equations are equivalent, it leaves 
open the question as to why this is so. By re-deriving the equations of motion using 
symplectic reduction, not only do we obtain both sets of equations, but the map S 
also follows naturally. 



3 The Neumann Connection 

The bundle vr : Q — > SE(2) is equipped with a principal fibre bundle connection, 
called the Neumann connection by Montgomery [1986]. There are many ways of 
describing this connection, but from a physical point of view, the definition using 
the horizontal lift operator is perhaps most appealing. From this point of view, the 
Neumann connection is a map 

h : Q X rSE(2) ^ TQ, h{ip, g, g) = o ip, (3.1) 

where $ is the solution of the Neumann problem (2.16) associated to {g,g). In other 
words, the Neumann connection associates to each motion {g,g) the corresponding 
induced velocity field of the fluid, and hence encodes the effect of the body on the 
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fluid. It is important to note that the Neumann connection does not depend on the 
point vortex model and is vahd for any vorticity field. 

Similar connections as this one have been described before (see for example Lewis, Marsden, Montgomery, ai 
[1986]; Khesin and Lee [2008]) but a complete overview of its definition and proper- 
ties seems to be lacking. In this section, we give an outline of the properties of the 
Neumann connection which are relevant for the developments in this paper, leaving 
detailed proofs for the appendix. 

Invariance of the Kirchhoff Decomposition. Before introducing the Neumann 
connection, we prove that the velocity potential ^{g,g;x.) is left SE(2)-invariant, 
expressing the fact that the dynamics is invariant under translations and rotations 
of the combined solid-fluid system. 

Proposition 3.1. The velocity potential <I> is left SFi{2) -invariant in the sense that 

<^{hg,TLh{g);hx) = ^{g,g;x) (3.2) 
for all h G SE(2) and (^,5;x) G rSE(2) x M^, 

Proof. This assertion can be proved in a number of different ways. The easiest 
is to use the assumption that the body is circular and solve the equation for the 
elementary potentials explicitly. Recall that these elementary potentials are given by 
(2.18). It is then straightforward to check that (3.2) holds, using the transformation 
properties of the velocity in the inertial frame. ■ 

The Connection One-form. For our purposes, it is convenient to define the 
Neumann connection through its connection one-form A : TQ jtvoi given by 

■^{ip,g){9,^) = ^*Uv, (3.3) 

where Uv is the divergence-free part in the Helmholtz-Hodge decomposition (2.15) 
of the Eulerian velocity u = 99 o cp^^. A proof that this prescription determines a 
well-defined connection form can be found in the appendix, proposition A.l, where 
it is also shown that this prescription agrees with the horizontal lift operator (3.1). 

The Curvature of the Neumann Connection. It will be convenient in what 
follows to have an expression for the ^u-component of the curvature of the Neumann 
connection, where for now fi is an arbitrary element of 3i*^^. Later on, /i will be the 
vorticity (2.20) associated to N point vortices. 

The curvature of a principal fiber bundle connection is a two-form B whose definition 
is listed in the appendix. For the Neumann connection we have 

%^)((5i,9ii),(52,9i2)) = -A<„^)([X^,y^]), (3.4) 

where and are horizontal vector fields on Q such that X^(g, 99) = {gi,ipi) 
and Y^{g,ip) = {g2,^2)- 
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Proposition 3.2. Let {gi,(pi) and {g2,'f2) be elements ofT(^g^^-jQ and denote the 
solutions of the Neumann problem (2.16) associated to {g,gi) resp. {g,g2) by 
and ^2- Then the ^-component of the curvature B is given by 

(^,S(<,,^)((5i,93i),(52,9^'2))> = ((Ai,d$i Ad$2))^- / a A*(d$i Ad$2), (3.5) 

where {{■,■)) is the metric on the space of forms on !F , induced by the Euclidian 
metric on T . 

Proof. Let be equal to (/?*da, and pick a such that a = u^, where Uv is a 
divergence-free vector field on T tangent to dT . For example, in the case of the 
vorticity due to N point vortices, Uv is the velocity field due to the vortices and their 
images (see Saffman [1992]). The calculation of the curvature involves computing 
the Jacobi-Lie bracket of two horizontal vector fields and taking the divergence-free 
part of the result. Because of the special form of a we can dispense with the latter 
step, since a is chosen to be L2-orthogonal to gradient vector fields. Therefore, the 
//-component of the curvature is given by 

(/i,S(,,^)((9l,^l),(52,9i2))> = -(/i,^(g,^)([X^,l^^])> 

(9.*a)([X^,y^])%, 



where the bracket on the left-hand side is the Jacobi-Lie bracket, which is the 
negative of the usual commutator of vector fields. 

The bracket can be made more explicit by noting that (as vector fields on Q) 

[x^, y^](9, V.) = ([(d<^>l)^ (da>2) V v', • • • ), 

where the dots denote a term in Tg SE(2) whose explicit form doesn't matter. The 
curvature then becomes 

(/., 0(3,^)((5i, 9^1), (52, <^2))> = - ^ a([(d$l)^ (dcl>2)1) dx. 

The remainder of the proof relies on the following formula for the co differential of 
a wedge product (Vaisman [1994, formula 1.34]): 

5(a A /?) = 5a A /? + {-l^a A 5/3 - [a^ /J^ , (3.6) 

where p = dega. Applying (3.6) with a = d<l>i and /3 = d<I>2 gives 

[(d^>l)^ (d$2)'']* = A$id$2 - A$2d$i - (5(d$i A d$2), 

where A is the Laplace-Beltrami operator on differential forms. Since both $i and 
$2 are harmonic, the first two terms of the right-hand side vanish. Consequently, 
the /i-component of the curvature becomes 

(/^,%,v)((5i,'/ii),(ff2,vi2))> = ((a,5(d$i Ad$2)))^ 

= ((da, d<I>i A d$2))jr - / a A *(d$i A d$2), 

JdT 
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using the adjointness of d and 5 for manifolds with boundary (this is a simple 
consequence of Stokes' theorem). ■ 

In traditional fluid mechanics notation, formula (3.5) becomes 

{li,B) = I (V X uy) • (V$i X V^2)dx - I uy • (n X (V$i X \I^2))dA, 
a formula first derived in Montgomery [1986]. 

4 Reduction with Respect to the Diffeomorphism Group 

As mentioned in the introduction, the fluid-solid system on the space Q C SE(2) x 
Emb(^05]^^) is invariant under the particle relabelling group Diffvoi- We now per- 
form symplectic reduction to eliminate that symmetry. In order to do so, we need 
to fix a value of the momentum map associated to the Diffyoi-symmetry. From a 
physical point of view, this boils down to fixing the vorticity of the system; it is at 
this point that the assumption is used that the vorticity is concentrated in point 
vortices. 

Before tackling symplectic reduction in the context of the fluid-solid system, we first 

give a general overview of cotangent bundle reduction following Marsden, Misiolek, Ortega, Perlmutter, and 

[2007]. Roughly speaking, applying symplectic reduction to a cotangent bundle 

yields a space which is diffeomorphic to the product of a reduced cotangent bundle 

and a co-adjoint orbit of the group. The dynamics on the reduced space is governed 

by a reduced Hamiltonian and a magnetic symplectic form: the reduced symplectic 

form is the sum of the canonical symplectic forms on the individual factors and an 

additional magnetic term. 

As we shall see below, in the case of the fluid-solid system, the reduced phase space 
is equal to T* SE(2) x M^^, where the first factor describes the rigid body, while 
the second factor determines the configuration of the vortices. At first sight, it may 
therefore appear that the intermediate fluid is completely gone. Yet, the vortices 
act on the rigid body, and vice versa, through the surrounding fluid. The answer 
to this apparent contradiction is that the effect of the fluid is concentrated in the 
magnetic symplectic form, for which we derive a convenient expression below. 

4.1 Cotangent Bundle Reduction: Review 

In this section, we collect some relevant results from Marsden, Misiolek, Ortega, Perlmutter, and Ratiu 

[2007]. We consider a manifold Q on which a Lie group G, with Lie algebra g, acts 

from the right and we denote the action by cr : Q x G —> Q. In addition, we assume 

that we are given a connection one-form A : TQ q with curvature two-form B. 

In the rest of this paper, Q will be the configuration space of the solid-fluid system, 

while the structure group G will be Diffvoi, and the connection A the Neumann 

connection. 
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The Curvature as a Two-form on the Reduced Space. Let /U be an element 
of Q*, and denote its isotropy subgroup under the co-adjoint action by G^: 

Gf, = {g£G: CoAdgfi = ^i} . 

Consider the contraction {fj,, B) of with the curvature B: due to the G-equivariance 
and the fact that B vanishes on vertical vectors (see (A. 2)), we may show that (/u, B) 
is a G^-invariant form. 

Proposition 4.1. The ^-component of the curvature B is a G ^-invariant two-form: 

1. a* {fi, B) = {fi, B), for all g G G^; 

2. ii-Q {fi,B) = for all^e Q^,. 

Hence, {fJ^,B) induces a two-form B^ on Q/G^ such that 

^*Q,G^^^^ = il^^lS), 
where t^q,g^_, ■ Q Q/G^ is the quotient map. 
Proof. For the first property, we have 

a; {f,,B) = {i^,a;B) = {fi,Adg-^B) = (CoAdg(^),e) = {is,B) 

if 5 G G^. The second item follows from the corresponding property for B. Actually, 
more is true: i^g B) = for all ^ £ Q, not just in g^. ■ 

Cotangent Bundle Reduction. The first stage in reducing the phase space con- 
sists of dividing out the particle relabelling symmetry. To this end, we use the frame- 
work of cotangent bundle reduction to construct the reduced phase space. The 
framework for cotangent bundle reduction outlined in Marsden, Misiolek, Ortega, Perlmutter, and Ratiu 
[2007] allows us to write down the reduced phase space and the modified symplectic 
form. We quote: 

Theorem 4.2. (Marsden et al. [theorem 2.2.1]) 

1. There exists a symplectic imbedding of the reduced phase space {T*Q)^ into 
the cotangent bundle T*{Q/G^) with the shifted symplectic structure := 

^can . 

2. The image of ip^ is the vector subbundle [TTTQ^Gfj,iV)]° ofT*{Q/G^), where 
V C TQ is the vector subbundle consisting of vectors tangent to the G-orbits 
in Q, and ° denotes the annihilator relative to the natural duality pairing 
between T{Q/G^) andT*{Q/G^). 
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Here {T*Q)^ is the reduced phase space J^^{n)/G^, where is the isotropy sub- 
group of II. The two- form is called the magnetic two-form^ and is defined 
as = T^QiQ (3^, where t^q/g^ '■ T*{Q/G^) — > Q/G^ is the cotangent bundle 
projection and is a two- form on Q/Gfj,. 

To define (3^, consider the one-form := on Q. One can show that da^ 

(but not itself) is G^-invariant, and the induced two-form on Q/G^ is precisely 
or in other words 

where ttq^g^ : Q Q/G^ is the quotient map. The calculation of da^ is greatly 
facilitated by using the Cartan structure equation (Kobayashi and Nomizu [1963]), 
which allows us to rewrite da^ as 

da^ = {^x,B)-{^l,[A,A)■ (4.1) 

Observe the sign difference on the right-hand side with Marsden, Misiolek, Ortega, Perlmutter, and Ratiu 
[2007], which is due to the fact that we take G to be acting on Q from the right. 

4.2 The Reduced Phase Space 

Theorem 4.2 provides us with an explicit prescription for the reduced phase space. 
Recall that the diffeomorphism group Diffvoi acts on the dual Lie algebra dr2^(.7^o) 
through pull-back. In particular, if /i represents the vorticity due to N point vortices 
as in (2.20), then 

N 

(t>*n = Ti6{5i - (t>-\5ii)) dx A dy. 

i=l 

In other words, the diffeomorphism group acts on the space of point vortices by 
simply moving the vortices. It also follows that the isotropy subgroup of /i consists of 
all diffeomorphisms (j) for which the vortex reference locations are fixed: 0(xi) = Xi, 
for i = 1, . . . , N . We denote the group of all such diffeomorphisms as Diff^oi,^. 

The group Diffvoi,/^ acts on Q and moreover, the quotient of Q by this action 
is diffeomorphic to SE(2) x M^^. To see this, note that Diffvoi,/^ acts on the 
Emb(.Fo,M^) factor only, and that there exists a diffeomorphism of the quotient 
space Emb(j^o,IK2)/Diff™i,^, with M^JV, given by 

[if] ^ {ip{5ii),...,(f{xN)). 

Similarly, the projection of Q onto the quotient space Q/Diff voi,/i = SE(2) x M.^^ is 
given by 

^Q,G^ ■■ {g, ^ (a; fi^i), 'p{^n))- 

Proposition 4.3. After reducing by the group of volume preserving diffeomor- 
phisms, the reduced phase is given by T* SE(2) x M^^. 
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Proof. The vertical bundle V on Q consists of vectors {g,0;(p,<f). Projecting this 
bundle down under Tttq^g^ shows that Tttq^g^{V) is spanned by elements of the 
form 

T-n-Q,G^{g,0;ip,(p) = (5',O;v?(xi),0(xi); . . . ; (^(xat), (^(xat)) 

and as each of the N vectors on the right hand side can range over the whole of TM^, 
T7TQ^Gf,{V) is equal to SE(2) xTM^A'. Its annihilator is therefore T* SE(2) xM^A'. ■ 

The reduced symplectic structure on T* SE(2) x M?^ is described in theorem 4.2. 
Explicitly, it is given by 

= ^can — B^, (4.2) 

where ilcan is the canonical symplectic structure on T* SE(2) and is the pullback 
to T* SE(2) X M^Af of the form (3^ on SE(2) x M^a^. Prom now on, we will no longer 
make any notational distinction between (3^ and its pull-back and denote both 

Apart from the reduced symplectic form, which will be determined explicitly later 
on, the dynamics on the reduced space is also governed by a reduced version of the 
Hamiltonian (2.13). The calculation of this Hamiltonian is the subject of the next 
section. 



4.3 The Reduced Hamiltonian 

The kinetic energy (2.6) is invariant under the action of Diffvoi,^ on Q and hence 
determines a reduced kinetic energy function on the quotient space SE(2) x M^^. 
Parts of the computation of the explicit form of the reduced kinetic energy can be 
found throughout the literature, but since no single reference has a complete picture, 
we briefly recall these results here. 

Using the Helmholtz- Hodge decomposition (2.15), the kinetic energy (2.6) of the 
combined solid-fluid system can be written as 

1 If 

T{9, "P) = 2 5'))sE(2) + 2 u • u dx 

= I {{9, 5))sE(2) + ^ ^ uy uy dx + ^ ^ V$ ■ V$ dx, (4.3) 

where we have used the fact that uy and are L2-orthogonal. 

Following a reasoning similar to Marsden and Weinstein [1983], it can be shown 
that the kinetic energy of the vortex system is nothing but the negative of the 
Kirchhoff-Routh function Wg (see (2.22)): 

uy • uy dx = -Wg(Xi, . . . 

The gradient term in (4.3) can be rewritten using a standard procedure, going back 
to Kirchhoff and Lamb, and yields the familiar added masses and moments of inertia 
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for a rigid body in potential flow. By using Green's theorem to rewrite the gradi- 
ent term as an integral over the boundary, and substituting the Kirchhoff expansion 
(2.17), one can show that (see Lamb [1945]; Kanso, Marsden, Rowley, and Melli-Huber 
[2005]) 

^^V$-V<I>dx = i(J7 V)M, (^^), (4.4) 

where is a 3-by-3 matrix of added masses and inertia whose entries depend only 
on the geometry of the rigid body. For the circular cylinder, IWq^ can be evaluated 
explicitly: 

The important point is to note that the gradient term (4.4) has the same form as 
the kinetic energy of the rigid body. By introducing the matrix MI — Mll^ -t- IVEjj, 
where is the mass matrix (2.7), the kinetic energy of the rigid body, together 
with the gradient term, can be written as 

which determines by left extension a function, also denoted by T, on TSE(2). 
Putting everything together, we conclude that the total reduced kinetic energy on 
T SE(2) X M^Af is given by 

Tred(f^,V;Xi,...,X^) = ^(17 V)m(^^^ -VFg(Xi,...,Xjv), (4.5) 

where Wg is given by (2.22) and (2.23). 

Two issues are noteworthy here. First of all, we follow Marsden and Weinstein 
[1983] and Shashikanth, Marsden, Burdick, and Kelly [2002] in "regularizing" Wg 
by subtracting infinite contributions that arise when putting Xq = Xi in the ex- 
pression for the Green's function G in (2.23). Secondly, we have chosen to express 
the kinetic energy in the body frame, which will be more convenient later on. 

4.4 The Magnetic Two-form 

Recall that is a two-form on X SE(2) defined as 

/3^(xi, . . . ,XAr;5r)((vi, . . . , VAr,s-i), (wi, . . .,^^,92)) = da^{ip, g){{ipi, gi), ((^2,52)), 

where Vj,Wj € Tx.M^ (i = 1, . . . ,N), and gi and g2 are elements of Tg SE{2). The 
embedding (p has to satisfy ¥'(xi) = Xj for i = 1,. . . ,N, while cpi and (^2 should 
satisfy 

V9i(xj)=Vj and (p2i5ci) = Wi (for i = 1, . . . , iV). 
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Using Cartan's structure equation (4.1), da^ can be rewritten as the difference of a 
curvature term and a "Lie-Poisson" term: if we introduce two-forms /3curv and (3\,p 
on M^Af X SE(2), determined by 

^Q,G^/5curv = (m, B) , and vr^,G^A.p = (/^, [-4, -^]) , 
then (3^ = Pcuw - A.p- 

In the computation of these terms, we wih frequently encounter expressions involving 
the vertical part of {ipi,gi) and {(p2,g2), evaluated at the vortex locations. We now 
derive a convenient expression for these terms. 

Let $1 be the solution of the Neumann problem (2.16) with boundary data gi and 
denote by u the Eulerian velocity field associated to (pi. From the Helmholtz-Hodge 
decomposition u = Uv + V<I>i, it follows the divergence-free part Uv satisfies (for 
i = l,...,N) 

Uv(Xi) = u(Xi) - V$l(Xi) 

= Vi - V«>i(xi), 

since u(xj) = ipi{5ii) = Vj. From the definition of A it then follows that we have 
proved the following fact about the vertical part of ipi: 

= r<^"i(vi- V$i(xi)) (4.6) 

for i = 1, . . . , N . A similar property holds for ip2, with Vj replaced by Wj and 
involving $2j the solution of (2.16) with boundary data g2- 

As an aside, we note that the Neumann connection induces a connection on the 
trivial bundle x SE(2) SE(2), and the expression between brackets on the 
right-hand side of (4.6) is the vertical component of the vector Vj. This connection 
is important in the theory of Routh reduction (see Marsden, Ratiu, and Scheurle 
[2000]). 

The Curvature Term. The curvature term /3curv is given by 

/?curv(xi,. . . ,XAr;S-)((vi,. . . ,V7V,5-i), (wi, . . .,WN,g2)) 

= {fJ',i3{(p,g){{^pi,gi), (^32, 52))) , 

where B is the curvature of the Neumann connection. It follows from (3.5) that 
PcuTv does not depend on the value of Vj and Wj, but only on gi and g2- From now 
on, we will therefore suppress these arguments. 

Proposition 4.4. For the rigid body interacting with N point vortices, the curvature 
term is given by 

N 

Pcurv (xi , . . . , X^v ; 9)(9i,92) = 5^ri(ia;(V$i(xi), V$2(x.)), (4.7) 

i=l 
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where $i and ^2 o,i"g the solutions of (2.16) associated to gi and §2, respectively, 
and dx is the volume element on T . 

Proof. We use the explicit form for the elementary velocity potentials derived in 
proposition 3.1. Note that /Jcurv depends on g\ and g2 through g~^g\ and g~^g2, so 
we can calculate its expression on the basis {e^,ex,ey} of 5e(2) given in (2.2) by 
assuming that g = e and 51, ^2 S SE(2) = se(2). 

First, we prove that the boundary term in (3.5) is zero. This is obvious if either gi 
or g2 is proportional to e^, since = 0. We may therefore assume that gi = 
and g2 = Gy. In this case, the two- form d^x A restricted to the boundary of 
(the curve (x — xq)^ + {u — Vof' = 1) is given by d^x A d^y = dx. Since *(ix = 1 in 
two dimensions, the boundary term reduces to 

/ a = r, (4.8) 
Jar 

which is zero, since by assumption there is no circulation around the rigid body. 

This leaves us with the integral over the fluid domain in (3.5). Plugging in the 
explicit form (2.20) of the vorticity, we have 

/3curv(xi, . . . ,yLN;g){gi,92) = ((^,d^>l A d^2))jr 

N 

= ^r,(ix(v^>i(x,),v$2(x^)). 

i=l 

We do not work out this expression any further, since it will be cancelled out by a 
similar term in the computation of /3i.p. ■ 

It is remarkable that the contribution from the curvature to the magnetic term 
will be cancelled in its entirety by a similar term in the expression for Lie-Poisson 
term, which is computed below. However, in the case of non-zero circulation, the 
curvature generates an additional term (4.8) proportional to the circulation. The 
effect of this term is studied in Vankerschaver, Kanso, and Marsden [2008]. 

The Lie-Poisson Term. The Lie-Poisson term /3i p is given by 

A.p(xi, . . . ,XAr;5f)((vi, . . . , VAr,gi), (wi,. . . , WAr,g2)) 

= [■A{v,g){^l, ai), ^{ip,g){^2, 92)]) , 

where {{pi^gi) and {^2^92) have similar interpretations as before. The right-hand 
side can be made more explicit by noting that, for divergence-free vector fields ui 
and U2 tangent to dJ- and arbitrary one- forms a, the following holds: 



/ a([ui, U2]) dx = / da(ui, U2) dx. 



(4.9) 
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The proof of this assertion is a straightforward apphcation of Cartan's magic formula 
and parallels the proof of theorem 4.2 in Marsden and Weinstein [1983]. 

Using this formula, we have 

N 

= ^^i {V*dx) {{A(^^^g){(pi,gi)){5ii), {A(^^^g){ip2,92)){^i)) ■ 
1=1 

By substituting (4.6) in this expression, we conclude that the Lie-Poisson term is 
given by 

A.p(xi, . . . ,XAr;5)((vi, . . .,VN,gi), (Wi,. . . , W7V,52)) 

N 

= ^ r, £ix(v, - V^>l(x,), W, - v$2(x0). 

i=l 

The Magnetic Two-form: Putting Everything Together. Using the previ- 
ously derived results, we may conclude that is given by (suppressing its arguments 
for the sake of clarity) 

— Pcurv Pl.p 
N 

= ri(dx(V^>l(Xi), V$2(Xi)) - dx(v, - V$l(Xi), Wi - V$2(Xi)) 

i=l 

N 

= Y^i{- dx(vi,Wi) + dx(vi, V^'2(xi)) +(ix(V^>i(xi),Wi)). 

i=l 

Introducing stream functions ^' ■(5, g; x) (i = 1, 2) as harmonic conjugates to ^i{g, g; x), 
i.e. such that 

V$Jdx = d^'i i = l,2 (4.10) 
we may rewrite the expression for f3^ as 

N 

f3^ = Y,ri{- dx(vi, w,) - d^^(v,) + d^;(w,)) . 

i=l 

The stream function can be written in Lamb form by introducing the elementary 
stream functions ^a;, and as harmonic conjugates to $2;, ^y, and re- 
spectively. For i = 1,2, we then have 

'^'1(9,9;^) = Ux,i*x(5',x) -h%i^j^(5f,x) + Va;,j^cu(5,x) 

= dxo{gi)^x{g,^) +dyo{gi)'^y{g,x) + d(3{gi)'^^{g,x), 
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and so the magnetic two-form becomes 

N 

1=1 

where the coordinates on SE(2) are denoted by = {xo,yo, 13}, and ^'^ = {^^, ^' 

Theorem 4.5. The magnetic two-form (3^ is a two-form on M^^ X SE(2) and can 
be written as 

/3^(xi, . . . ,XAr;5)((vi, . . . ,VAr,5(i), (wi, . . .,^^,92)) 

N 



= (-da;(vi,Wi) + de(5, Xi)(5i, v^; 52, w^)) . (4.11) 

j=i 

where Q is the one-form on x SE(2) given by 

e(<7,x)(g,v) = ^'Ag,^)dU{9). (4.12) 

Proof. Note that the exterior derivative on the right-hand side of (4.11) is the 
exterior derivative on X SE(2). Expanding the second term, we have 

d{^'A{9,^)dU) = d^'^AdU 

= d^'^ A dxo -\- d^y A dyo 

= d]K2^'^ A dxo + djs^2'^y A dyo + ( ^ - ^1 dxo A dyo, 

where the subscript 'M^' indicates that the derivative should be taken with respect 
to the spatial variables only, i.e. that the SE(2)-variable should be kept fixed. Also, 
we rely on the fact that ^'^ = (since <I>(^ = 0) , and that ^'^ y does not depend on 
p. 

The theorem is proved once we show that the term in brackets vanishes. To show 
this, we need the explicit form for the elementary stream functions: 

^' = y-lJ^ and = (4 13) 

^ (x - Xo)2 + (y - yo)2 y - xqY + {y - yoY 

and it is then easily shown that this is indeed the case. ■ 

Putting everything together, we conclude that the symplectic structure Q^is on m?^ X 
T* SE(2) (see (4.2)) is given by 

= ^c.n - f3^,, (4.14) 

where is given in the theorem above. The first term is the canonical symplectic 
form on T* SE(2) while the magnetic consists of two parts, one (proportional to 
(ix) being the symplectic form on M."^^ , while the other one encodes the effects of 
the ambient fluid through the stream functions and "ify. Note that is a 

CO- adjoint orbit of the diffeomorphism group, and that its symplectic structure is 
nothing but the Konstant-Kirillov-Souriau form. 
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5 Reduction with Respect to SE{2) 

After reducing the fluid-solid system with respect to the Diffvoi-symmetry, we are 
left with a system on T* SE(2) x M?^ with a magnetic symplectic form. The group 
SE(2) acts diagonally on the reduced phase space. What remains is to carry out 
the reduction with respect to the group SE(2) and then to make the connection 
with the equations of motion in Shashikanth, Marsden, Burdick, and Kelly [2002] 
and Borisov and Mamaev [2003]. 

Two factors complicate this strategy, however: first of all, there is the presence of the 
magnetic term in the symplectic structure, which precludes using standard Euler- 
Poincare theory, for instance. Secondly, the phase space is a Cartesian product on 
which the symmetry group acts diagonally. Even in the absence of magnetic terms, 
this would lead to the appearance of certain interaction terms in the Poisson 
structure (see Krishnaprasad and Marsden [1987]). 

In order to do Poisson reduction for this kind of manifold, we begin this sec- 
tion by considering the product of a cotangent bundle T*H, where H is a Lie 
group, and a symplectic manifold P, and we assume that the symplectic struc- 
ture on T*H X P is the sum of the canonical symplectic structures on both fac- 
tors plus an additional magnetic term. This case extends both magnetic Lie- 
Poisson reduction in Marsden, Misiolek, Ortega, Perlmutter, and Ratiu [2007], for 
which the P-factor is absent, as well as the "coupling to a Lie group" scenario of 
Krishnaprasad and Marsden [1987], for which there is no magnetic term. Since the 
proofs in this section are rather lengthy, we have relegated them to Appendix B and 
we simply quote the relevant expressions here. 

The bulk of this section is then devoted to making these results explicit for the case 
where H = SE(2), P = M?'^ and the symplectic form is the magnetic symplectic 
form derived in the previous section; see equations (4.11) and (4.14). The induced 
Poisson structure with interaction terms will turn out to be nothing but the BMR 
Poisson structure (2.30). Secondly, the Poisson map induced by the momentum map 
J is just the Shashikanth map described in Theorem 2.2, so that the transformed 
Poisson structure is the SMBK Poisson structure (2.27). 

5.1 Coupling to a Lie Group 

As mentioned above, the purpose of this section is to generalize the reduction of 
phase spaces of the form T*H x P, where H is a Lie group acting on the symplectic 
manifold P. This is more general than results in the literature, since the symplectic 
structure onT*HxP may contain magnetic terms. We shall need this generalization 
in our application. The proofs of the theorems in this section can be found in 
appendix B. 

Specifically, if ujp denotes the symplectic form on P, we assume that the symplectic 
structure Q on T*H x P is given by 



(5.1) 
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where Wcan is the canonical symplectic structure on T*H and B is a closed two- form 
on G X P such that B^p = —ujp, i.e. the restriction of to P is just — cjp. The 
latter ensures that O is a well-defined symplectic form on T*H x P. 

We assume furthermore that H acts on P and we denote the action oi h £ H on 
p £ P hy h ■ p. Furthermore, we assume that this action leaves invariant top. The 
group H then acts diagonally on T*H x P: \i h £ H and (ag,p) G T*H x P, then 

h-{ag,p) = {T*Lf,-i{ag),h-p). (5.2) 

The Momentum Map. Consider an element ^ of f). We denote the fundamental 
vector field associated to on r*if x P by 

^T''Hxp{ag,p) := ^exp(tO • {ag,p) 

and similar for the fundamental vector field ^hxP on H x P. Under certain topo- 
logical assumptions (which are satisfied for the solid-fiuid system), there exists a 
momentum map J : T*H x P ^ I)* for the action of H, defined by 

where = (J, 0- It follows from the explicit form for the symplectic form (5.1) 
that J can be written as J = Jean — where Jean : T*H — s- ()* is the momentum 
map associated to the canonical symplectic form on T*H: 

igj^^Wcan = d (Jean, (5-3) 

while (j) : H X P ^ [)* is the so-called BQ-potential (borrowing the terminology of 
Marsden, Misiolek, Ortega, Perlmutter, and Ratiu [2007]): 

k„^,B = d{c^,0. (5.4) 

The theory of ;Bg-potentials and Poisson reduction for cotangent bundles T*H en- 
dowed with a magnetic symplectic form is further developed in Marsden, Misiolek, Ortega, Perlmutter, and 
[2007]. Here, we are dealing with a product T*H x P, for which no such theory can 
be found in the current literature. 

Infinitesimal Equivariance of the Momentum Map. Whereas the canonical 
momentum map Jean is equivariant, the same does not necessarily hold for the 
;B0-potential. To measure non-equivariance, we introduce a one-cocycle cj : iJ — > 
L{^,C{H X P)) by defining first a family of functions F^^^ (where r] £\] and g G H): 

'^r,,g{Kp) = - {(l){gh,gp),r]) + (Ad*(/)(/i,p), r/) , 

for all {h,p) £ H X P, and then putting cr{g) ■ ij = F^^^. This definition follows the 

usual introduction of cocycles for momentum maps; see for instance Guillemin and Sternberg 

[1984]; Marsden and Ratiu [1994]. 
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The one-cocycle a induces a two-cocycle S : f) x f) — > C'^{H x P) given by 

where an '■ H ^ C°°{H x P) is defined by Unig) = (j{g) ■ rj. It is not hard to verify 
that Ti is exphcitly given by 

r]) = - (0, 7?]) + Bi^HxP, mxp). (5.5) 



The Poisson Structure on i)* x P. The canonical Poisson structure on T*H x P 
associated to the symplectic structure (5.1) gives rise to a Poisson structure on the 
quotient space f)* x P, which we denote by {•, - jint- 

The exphcit form of this Poisson structure is derived in Appendix B. Before quoting 
this result, we first introduce an operation * : x f) — > C°°{H x P) defined by 

where {■,-}p is the Poisson structure associated to ujp, and an operation < : x 
C°°{P) ^C^iH X P) by putting 

e<]F:={</.^|P,F}p. 

Using these two operations, the reduced Poisson structure is given in the following 
theorem. Note the interaction term due to curvature and the last two terms which 
are due to the coupling of the Lie group H with the symplectic manifold P. 

Theorem 5.1. The reduced Poisson structure on i)* x P is given by 
for functions f = f{fi,x) and k = k{i^i,x) on f)* x P. 



Shifting Away the Interaction Terms. The reduced Poisson structure {•, - jint 
described in theorem 5.1 is not canonical, i.e. it is not the sum of the Lie- Poisson 
structure on f)* and the canonical Poisson structure on M^'^, but rather contains a 
number of interaction terms. However, using the ,B0-potential 4>, we can define 
a shift map S from f)* x to itself, taking {•, - jint into the canonical Poisson 
structure. The price we have to pay for getting rid of interaction terms is the 
introduction of the non-equivariance cocycle S of (j), as in the following definition. 

Definition 5.2. The natural Poisson structure on t)* x P is given by 

= A:|^4i.p - {f\P,k\p}p - S ^) (5.7) 
for functions f = f{p,x) and k = k{fi,x) on t)* x P. 
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The main result of this section is then given in the following theorem. It related 
the Poisson structure with interaction terms with the natural one (possibly with a 
cocycle). The former will turn out the BMR Poisson structure, while the latter is 
nothing but the SMBK structure. The cocycle will encode the effects of nonzero 
circulation on the rigid body. 

Theorem 5.3. The map 5 : f)* x P ^ {)* x P given by 

S{fi,x) = {fi - (p{e,x),x). (5.8) 

is a Poisson isomorphism taking the Poisson structure {•, -jint with interaction terms 
into the Poisson structure 

{foS,ko 5}int = {/, k}j: o S. (5.9) 

The Symplectic Leaves. The shifted Poisson structure {•, -js is the sum of the 
Lie- Poisson structure on the canonical Poisson structure on P, and a cocycle 
term. If S = 0, this allows us to write down a convenient expression for the sym- 
plectic leaves in f)* x P: 

Proposition 5.4. For S = 0, the symplectic leaves inl)* xP of the Poisson structure 
{■)'}s=o fl'^e of the form x P, where is the co-adjoint orbit of an element 

The proof follows that of proposition 10.3.3 in Marsden, Misiolek, Ortega, Perlmutter, and Ratiu 

[2007] and relies on the fact that the symplectic leaves are precisely the symplectic 

reduced spaces. But since the Poisson structure on the reduced space is simply the 

sum of the Lie-Poisson and the canonical Poisson structure, the reduced space at /i 

is Of^ x P, as above. When S 7^ 0, we expect the co-adjoint orbit to be replaced 

by an orbit Oj^ of a suitable affine action of i/ on {)*. 

5.2 The Fluid-Solid System 

Now we are ready to specialize the theory in the previous sections to the case of 
a solid interacting dynamically with point vortices. Recall that our goal is to find 
the Hamiltonian structure for this problem using reduction techniques and to use 
the shifting map developed in the preceding section to relate the two Hamiltonian 
structures in the literature. 

Recall that the reduced phase space for the solid-fluid system is T* SE(2) x M^a^. 
Now, the group SE(2) acts on SE(2) x by the diagonal left action, denoted by 
$ and given in inertial coordinates by 

$/i(g;xi, . . . ,X7v) = {hg;hxi, . . . ,hxN)- 

Hence, SE(2) acts from the left on T* SE(2) x M^A^ using the cotangent lift in the 
first factor, and thus, this action leaves the Hamiltonian (4.5) invariant. 
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This system is of the form T*H x P, as discussed in the previous section, where 
H = SE(2) and P = M^^. We now apply the results of that section to divide out 
the SE(2)-symmetry and obtain a system on se(2)* x M."^^ . 

This reduction is similar to the passage from inertial to body coordinates for the 
rigid body (see for example Marsden and Ratiu [1994]). To see this, notice that the 
space 5e(2)* x M^^ is obtained from T* SE(2) x M^^ by dividing out by the diagonal 
SE(2)-action: the quotient mapping is given by 

p : T* SE(2) x M^^ ^ se(2)* x R^^ , 

defined as 

p{g,ag;xi,. . . ,xn) = {T*Lg{ag);Xi, . . . 
where, ii g = (ii, xq), then Xj and Xj are related by 

Xi = RXi + xo. (5.10) 

In other words, if Xj describes the location of the ith vortex in inertial coordinates, 
then Xj is its location in a frame fixed to the body. 

Proposition 5.5. The magnetic symplectic structure (4.14) is invariant under the 
action of SE(2) on T* SE(2) x M^^ described above. 



Proof. Recall the expression (4.14) for rjg, where is given by (4.11). It is 
a standard result (see for instance Cushman and Bates [1997]) that the canonical 
symplectic form ilcan on T* SE(2) is invariant under the left action of SE(2), and a 
similar result holds for the form "^^i rj/i(vj, Wj). The only thing that remains to 
be shown is that is SE(2)-invariant, but this follows from the SE(2)-invariance of 
^, which is itself a consequence of proposition 3.1 and the fact that ^ and ^ are 
harmonic conjugates. ■ 



The Momentum Map. Recall from section 5.1 that the momentum map J for 
the action of SE(2) on T* SE(2) x is the difference of two separate parts: J = 
Jean — where Jean is the momentum map (5.3) defined by the canonical symplectic 
form on T* SE(2), while (j) is the so-called i3g-potential (5.4). 

The momentum map J is a map from T* SE(2) x to se(2)*. It will be convenient 
to identify T* SE(2) with SE(2) x se(2)* through left translations, and to use the 
fact that se(2)* is isomorphic to M'^, so that J is a collection of three functions 
{J^,jy,J^) on M^^ X SE(2) x se{2)* . A typical element of that space will be 
denoted as (xi, . . . g,Il), where g = (i?, xq) and 11 = (11^., IIj^, IIj^), but for the 
sake of clarity we will usually suppress the argument of J. 

Proposition 5.6. The momentum map J associated to the SFi{2)-symmmetry repre- 
sents the spatial momentum of the solid-fluid system, and is given by J = (J^, Jy, Jq), 
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where 

(t) - « (":) - g -^-p' ^ - "> 

i=l 1=1 
where F = X^^^i Tj is i/ie toia/ vortex strength. 

Proof. The canonical part Jean can be obtained through a standard calculation for 
cotangent lifted actions (see Marsden and Ratiu [1994]). The result is 

for g = {R,xo) and 11 = (11^,11^,110). 

For the jBg-potential, note that the form is the sum of a "pure vortex" part {i.e. 
not involving the fluid) and a part involving the stream functions {i.e. the form B). 

The i3g-potential corresponding to the pure-vortex part is given by 

/AX \ ^ ^ T 

^vortex\ and cpl,,,^ = Y,^{^I + yf)- 

VPvortox/ .^^ ^ 

These expressions coincide up to sign with those in Adams and Ratiu [1988], which 
is a consequence of the fact that our symplectic structure is the negative of theirs. 

Finally, for the stream function term dQ in observe that G is an SE(2)-invariant 
one-form on SE(2) x M^^. The jBg-potential associated to d© is hence given by (see 
the remark following Theorem 7.1.1 in Marsden, Misiolek, Ortega, Perlmutter, and Ratiu 
[2007]) 

((^stream, = ksEi2)xR2N ® ^^"^ ^ ^ 

Explicitly, 

ULeaJ-|^Ef=ir.^,(,,X.)J ^st.ean.-l^r.^..(,,X.), 

where the elementary stream functions are given by (4.13). 

The full momentum map J is then the sum of these three contributions: 

J — Jean '/'vortex '/'stream j 

and this is precisely (5.11). ■ 
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The expression for the momentum map can be rewritten by introducing the mo- 
mentum map of the soUd-fluid system in body coordinates (Shashikanth [2005]; 
Kanso and Oskouei [2008]): 



N N 
i=l i=l 

where ^x,^Y,^n are the expression for the elementary stream functions in body 
coordinates. By using the relation (5.10) between inertial and body coordinates, 
one can show that the spatial and body momentum are related by 

=r(^^J^^ +T^oxe3 and = Jf, + (xq x J,) • eg, (5.12) 



where Jx = {Jx, Jy)^ ■ Note that the spatial momentum is a function on SE(2) xl 
whereas the body momentum is a function on M^^, and that the SE(2)-dependence 
of J is determined by the relation above. In particular, if we evaluate J at the iden- 
tity, putting R = 1 and xq = in (5.12), then we obtain just the body momentum: 

J(e,x) = J(X). 



Non-equivariance of the Momentum Map. Using the definition (5.5), the 
non-equivariance two-cocycle of the momentum map is a map T, : se(2) x se(2) — > 
C°°(SE(2) X M^^). To compute S, we recall the basis {en,e^.,ej^} of se(2) given in 
(2.2) and let {e^,e*.,e*} be the corresponding dual basis of se(2)*. 

Theorem 5.7. The non-equivariance two-cocycle S of the momentum map J is 
given by 

E = -re;Ae*, 
where F = Xl^^x total vortex strength. 

Proof. In order to use (5.5) to compute S, we need the infinitesimal generators 
corresponding to the basis elements e^,, e^ and e^. The infinitesimal generator of 
any element ^ G 5t{2), evaluated at (e; xi, . . . , x^r), will be denoted by ^ and is given 

by 

C '■= ?SE(2)xR2J^(e;Xi, . . . ,XAr) = (^; ^K2iv (xi) , . . . , ^ir2JV (xat)), 

where ^^2n is the infinitesimal generator of for the fundamental action of SE(2) 
on M^. Explicitly, 
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Evaluating on these vectors, we obtain 



-r 



/3/.(ex,en) = ( -X, + 

^ ( 1 




j=i 




Remarkably, when calculating S, the last two expressions are cancelled entirely by 
opposite contributions from the remaining terms: 



and similar for S(ey,en). The only non-zero term is given by Yj(<s.x^Gy) = — F. ■ 

In its current form, S has the same form as the non-equivariance two-cocycle for 
the A^- vortex problem in an unbounded fluid (see Adams and Ratiu [1988]). 

5.3 Poisson Structures 

Now we come to the conclusion of this paper. Using the theory developed in the 
previous sections, we derive an explicit form for the reduced Poisson bracket on 
se(2)* X associated to the magnetic symplectic form (4.14): this will turn out 
to be the BMR bracket. In the terminology of section 5.1, this is the bracket {•, - jint 
with interaction terms. Secondly, we then use the shift map S (see (5.8)) associated 
to the momentum map (5.11) to obtain the Poisson structure {•,•}£ where the 
interaction terms are absent, at the expense of a non-equivariance cocycle. When 
made more explicit, the latter bracket turns out to be the SMBK bracket. 

The BMR Poisson Structure. The bracket obtained by Poisson reduction of 
the magnetic symplectic structure on T* SE(2) X M^Af was described in theorem 5.1. 
The explicit computation of this Poisson bracket boils down to substituting the 
explicit expresssion (5.11) for the momentum map into the Poisson bracket {•, •} in 
theorem 5.1. After a long, but straightforward calculation, one obtains the BMR 
bracket (2.30). As an illustration, we compute one bracket element, leaving the 
others to the reader. 

We consider the functions F = Ha; and G = 11^ on se(2)* x and compute 
{IIj;, Ilylijit. Note that, considered as elements of se(2). 



S(ea;,en) 



- (0, [e^,en]) +/?^(ea;,en) 
0y + Pfj.{ex,en) 



= 




= e. 



■y 
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Computing (2.30) term by term, we have first of all that 

^ llx-ii'^-i?^ 

Tlx * n?/ — {4^xj '/*j/}vortex — ^ ^ ^ ) 

1=1 II Ml 

where {•, - jvortcx is the vortex bracket defined in (2.29). Secondly, 

(Gx ) ) = 

while the other terms are zero. Finally, it follows that 



N 

{n^,ny}int = r-^r, 

i=l 




which is precisely the first bracket element of (2.30). The computation of the other 
elements is similar. 

The SMBK Poisson Structure. By subjecting the BMR Poisson structure to 
the shift map S, we can eliminate the interaction terms from the Poisson structure, 
at the expense of introducing a non-equivariance cocycle. As shown in Theorem B.4, 
the result is a Poisson structure consisting of the sum of the Lie-Poisson and the 
vortex Poisson structure on the individual factors, together with a cocycle term: 

{F, G}s = {F|s5(2)* , G\g^(^2)* }l.p + {-^|IR2]V , G\^2N jvortex - S 

The last term is explicitly given by 

/5F 6F\ _ _^ r dF dG _ dG dF \ 

In the case where the total vorticity F is zero, this term vanishes and the bracket 
reduces to the SMBK bracket (B.4). 

The Symplectic Leaves. In the case where the two-cocycle S vanishes, or equiva- 
lently, when the total circulation F is zero, a convenient expression for the symplectic 
leaves in se(2)* x can be read off from proposition 5.4. 

Recall that the symplectic leaves for the Lie-Poisson structure in se(2)* = come 
in two varieties. One class consists of cylinders whose axis is the il-axis: 

O = {(n^., n.y, n) :nl + Ul = constant}, 

while the other class consists of the individual points (0, 0, O) of the O-axis. The 
symplectic leaves in se(2)* x are then the product of the symplectic leaves of 
the Lie-Poisson structure with M^^. 



/5F 6F\ 
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6 Conclusions and Future Work 

In this paper, we have used reduction theory to give a systematic derivation of the 
equations of motion for a circular rigid body in a perfect fluid with point vortices. 
Among other things, we have derived the Poisson structures that govern this prob- 
lem, and related them to the Poisson structures in the literature. However, the 
usefulness of our geometric method is not limited to merely reproducing the correct 
Poisson structure for one specific problem: by uncovering fundamental geometric 
structures such as the Neumann connection, we have shed new light on the precise 
nature of the solid-fluid interaction. Moreover, even though a number of simplifying 
assumptions were made in the beginning, it is expected that the present method can 
be extended without too much trouble to cover more general cases. Below, we have 
listed a number of open questions for which our method seems especially suited. 

Non-zero Circulation. As we have seen in this paper, if the total strength of 
the point vortices is non-zero, a cocycle term appears in the equations of motion. In 
a previous paper (see Vankerschaver, Kanso, and Marsden [2008]), we studied the 
case of a rigid body moving in a potential flow with circulation and found that the 
circulation manifests itself through the curvature of the Neumann connection. 

From a physical point of view, both phenomena are very similar. Indeed, having 
a non-zero circulation around the rigid body amounts to placing a vortex at the 
center of the rigid body. Hence, mathematically speaking, there appears to be a 
duality between the description in terms of cocycles and in terms of curvatures. 
Having a better mathematical understanding of this duality would prove to be use- 
ful for other physical theories as well, for example the dynamics of spin glasses in 
Cendra, Marsden, and Ratiu [2004]. In addition, it would be interesting to investi- 
gate the link between this theory with the work of Gay-Balmaz and Ratiu [2008]. 

Arbitrary Body Shapes. The pioneering work of SMBK was extended to cover 
the case of rigid bodies with arbitrary shapes by Shashikanth [2005], and studied fur- 
ther by Kanso and Oskouei [2008] with a view towards stability and bio-locomotion. 

By using methods of complex analysis, the methods of the present paper can easily 
be extended to cover the case of rigid bodies of arbitrary shape as well: any closed 
curve enclosing a simply-connected area of the plane can be mapped to the unit 
circle by a suitable conformal transformation and by pulling back the geometric 
objects in this paper along that map, arbitrary body shapes can be treated. 

Moreover, the objects defined in this paper transform naturally under conformal 
transformations. The magnetic form P^, for instance, will be mapped into a form 
which has the same appearance as before, but with now the elementary 
stream functions for a body of that shape. As known from the classical fluid dy- 
namics literature, these stream functions also transform naturally under conformal 
mappings. 
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Three-dimensional Bodies. In contrast to calculations using only standard, ad- 
hoc methods, our setup naturally generalizes to the case of three-dimensional flows 
interacting with rigid bodies. Shashikanth, Sheshmani, Kelly, and Marsden [2008] 
consider a rigid body interacting with vortex rings. The space of point vortices 
then is replaced by the space of vortex rings, and much of the analysis 
in the 2D case carries through, at the expense of significant computational work. 

From a geometric point of view, however, A4 plays a similar role as M^^: both are 
coadjoint orbits of the group Diffyoi of volume-preserving diffeomorphisms. Hence, 
Ai is equipped with a natural symplectic structure (see Arnold and Khesin [1998]), 
which is of importance for the calculation of (3^ in this case. Also, the definition 
of the Neumann connection remains valid in three dimensions. We therefore expect 
that the geometric approach can be extended without too much difficulty to this 
case, and will lead to a conceptually much clearer picture. 

Routhian Reduction by Stages and Dirac Reduction. It has been known 
for some time that the Lagrangian analogue of symplectic reduction is Routhian 
reduction (Marsden, Ratiu, and Scheurle [2000]). By starting from the Lagrangian 
of an incompressible fluid together with the rigid body Lagrangian, one could per- 
form reduction on the Lagrangian side using this theory. The advantage would be 
that Routhian reduction preserves the variational nature of the system. However, 
as the reduction procedure for the fluid-solid system consists of two succesive re- 
ductions, one would need a theory of Routhian reduction by stages. Developing 
such a theory would be of considerable interest. 

A related approach consists of using Dirac structures. Through the Hamilton- 
Pontryagin principle, Dirac structures offer a way of incorporating both the La- 
grangian and Hamiltonian formalisms making them ideally suited for systems with 
degenerate Lagrangians; see Yoshimura and Marsden [2006]. Since the point vortex 
Lagrangian is degenerate, Dirac structures are likely to be useful here. 

A Further Properties of the Neumann Connection 

Even though the Neumann connection is used in many problems from geometric 
fluid dynamics or differential geometry, a systematic presentation of its properties 
seems to be lacking. In this appendix we prove some basic theorems related to the 
Neumann connection. 

A.l Fiber Bundles and Connections. 

Recall that n : Q ^ Q/G is a principal fibre bundle with structure group G, then 
a connection one-form on Q is a g-valued one-form A on Q satisfying the following 
two properties: 
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1. A is G-equivariant: cr*A = Ad^-i o A, where ag : Q ^ Q denotes the 
G-action, and Ad is the adjoint action on g. 



2. Let be the infinitesimal generator associated to an element € g, i.e. 

d 

di' 



= ^o-(expt^,g) where a{g,q) = ag{q). 



Then A{^q) = t 

At each point q £ Q, the connection A induces complementary projection operators 
Pv,Ph -.TQ^TQ, given by 

Pvivg) = [A{vg)]Qiq), and Ph = 1-Pv, 

and referred to as the vertical and horizontal projection operators, respectively. 

More information on principal fiber bundles and connections can be found in Kobayashi and Nomizu 

[1963]. 



The Curvature. The curvature of a principal fiber bundle connection A is the 
0- valued two-form B defined as follows. Let Uq,Vq € TgQ and consider vector fields 
Xu and extending Ug and Vg, i.e. such that Xu{q) = Ug and Xy{q) = Vg. Denote 
the horizontal part of X^ as X^ := Pjj o Xu, and similar for X^ . Then, B can be 
conveniently expressed as 

Bg{Ug,Vg) = -Agi[X^,Y,'']). (A.l) 

With this definition, B is G-equivariant and vanishes on vertical vector fields: 

a*gB = Adg-iB and i^gB = 0, (A.2) 

for all (7 € G and G 0. 



A.2 The Neumann Connection 

The Connection One-form. We recall from section 3 that the connection one- 
form for the Neumann connection is given by 

where Uv is the divergence- free part in the Helmholtz-Hodge decomposition (2.15) 
of the Eulerian velocity u; see (3.3). 

Before showing that this prescription indeed yields a well-defined connection form, 
we note that each divergence-free vector field u G jC^oi defines a vertical vector field 
on Q, denoted by uq and given by 

"(3(5, 9^) = (0,T(^ o u). 

The vector field uq is the infinitesimal generator associated to u under the action 
of Diffvoi on Q. 
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Proposition A.l. The one-form A defined in (3.3) is a connection one-form. In 
other words, 

1. A is DiE^oi-equivariant: for all cp G Diffvoi o-'^T-d {g,g', ^, ^) S TQ, 

2. If u is an element of 3£voi with associated fundamental vector field uq, then 
A{uq) = u. 

Proof. To prove equivariance, note that the Eulerian velocity associated to o 
is equal to that associated to ip: 

u' = {(p o (p) o (^ip o <f>)~^ = ip o ip~^ = u. 

Hence, 

■^(g,<fio4>)ig, ipo(j)) = {Lpo 0)*Uv = </>*(99*Uv) = (l)*A(^g^^){g,ip). 

To prove the second property, note that the push- forward X = Lp^,\i \s divergence- 
free and tangent to the boundary of J- since u is divergence-free and tangent to dJ-Q 
and ip is volume-preserving. 

Evaluating the connection one-form on u therefore gives 

which concludes the proof that A is a well-defined connection one-form. ■ 

The Horizontal Lift. The horizontal lift operator associates to each {g, (p) a 
linear map h(g : Tg SE(2) — > T^^g^^)Q with the following properties: 

1. The composition T^^g^^^-K o h(g is the identity in Tg SE(2); 

2. The image of ii-{g.^) consists of horizontal vectors: A(^g^^) o h(g = 0. 

For the Neumann connection, the horizontal lift is given by 

h(9,^)(5) = (5,V$o99), 

where $ is the solution of the Neumann problem (2.16) associated to {g,g)] i.e. 
(w, v) = gg~^. It is easy to show that h satisfies both properties of the horizontal 
lift, and hence that h is indeed the horizontal lift of the Neumann connection as we 
claimed in section 3. 
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The Neumann Connection as a Mechanical Connection. When given a 
group-invariant metric on the total space of a principal fiber bundle, the mechanical 
connection is defined by declaring its horizontal subspaces to be orthogonal to the 
vertical bundle. The terminology stems from the fact that the underlying metric is 
usually the kinetic-energy metric of a mechanical system. 

Since there exists a natural metric (2.14) on Q, it should therefore come as no sur- 
prise that the Neumann connection can also be viewed as a mechanical connection. 

The horizontal subspace of the Neumann connection at an element (g, G Q is 
given by the kernel of Explicitly, a tangent vector {g,g; (p, (p) is in the kernel 

of if and only if ip = V<& o yj, where $ is the solution of the Neumann problem 

(2.16) associated to {g,g). 

On the other hand, the vertical subspace at {g, p) is generated by elements of the 
form VL.Q{g,p). Verifying that the Neumann connection is indeed the mechanical 
connection now boils down to checking that the vertical and horizontal subspaces 
are orthogonal with respect to the metric (2.14). 

Let (5, g] p, V<I> o p) and UQ(g, p) be horizontal and vertical tangent vectors, respec- 
tively. Then we have 



where the last equality follows from the L2-orthogonality between gradient vector 
fields and divergence-free vector fields tangent to the boundary of J^. 

B Proofs of the Theorems in Section 5.1 

In this technical appendix, we provide proofs for the theorems in section 5.1 regard- 
ing the reduced Poisson structures with or without interaction terms. We do so by 
using the diffeomorphism of T*H x P with H x t)* x P given by left translation. 
More specifically, we introduce a diffeomorphism A : T*H xP^Hx\j*xPhy 



and a diffeomorphism X: HxP—fHxP by X{g,p) = {g,g ^ ■ p) such that the 
following diagram commutes: 




V^>- ((/?*u)r/ = 0, 



M9,ag;p) = {g,T*Lg{ag),g ■ p) 



T*H X P 



A 



Hx\)* xP 



pr 



H xP 



X 



H xP 



Here, -kh ■ T*H x P H x P is given by Trnig, ag',p) = {g,p), while pr simply 
forgets the second factor: pv{g,fi,x) = {g,x). 
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The advantage of this construction is that under A the action (5.2) becomes left 
translation on the first factor. That is, if A{ag;p) = {g,fi,x), then A(/i • {ag;p)) = 
{hg,fi,x). The fact that the diagonal action (5.2) reduces in this way greatly sim- 
plifies Poisson reduction. 

B.l Symplectic Forms in Body Coordinates. 

We start by deriving in lemma B.l an explicit expression for the push-forward 

of the symplectic form. We then define a natural shift map S from H x i)* x P to 

itself, given by 

S{9, /U, x) = {g,^- (j){e, x),x), 

where : if x P — > ()* is the ;Bf|-potential associated to the action of H on HxP, i. e. 
the solution of (5.4). This map takes the symplectic structure (B.2) with interaction 
terms of lemma B.l into one where the interaction terms are absent. The explicit 
form of the resulting symplectic form is derived in lemma B.2. 

For future reference, we will refer to the push-forwards A^,u;can and A^fi as Wbody 
and r^bodyi respectively: 

t^body = A*CJcan and fibody = A^.^^. (B.l) 

Lemma B.l. The push-forward A^:Q of the symplectic form Q to H x i)* x P is 

given by 

{A^n){g,fi,x){{g ■ i,p,Vx), {g ■ ri,a,Wx)) = 

- B{e,x){S,Hxpie,x),r]Hxpie,x)) - B{e,x){{0,V:r), iO,Wx)) 

- (d05(e, x), (0, w^)) + {d(j)r,ie, x), (0, v,)) , (B.2) 

where {g,^,x) is an element of H x I)* x P and {g-S.^ p,Vx), {g-rj,a,Wx) are elements 
ofT^g,^,x){H X r X P) = TgH X r X TxP. 

Proof. Notice first that A^Q = A^cjcan — pr*(A*i?). The first term is the expression 
for the canonical symplectic form in body coordinates calculated by R. Cushman 
(quoted in Abraham and Marsden [1978], proposition 4.4.1) and is given by 

{A^uJciin){g,P',x){{g ■ ^,p,Vx), {g ■ ri,a,Wx)) = - {p,v) + {(^,0 + (/^> [^,v]) ■ 

For the second term, we have 

{KB){g,x){{g ■ ^,Vx), {g ■ v,Wx)) 

= B{X-\g,x)){TX~\g ■ ^,Vx),TX-\g ■ r,,Wx)) 

= Big,g- x){ig ■ C,T^givx) + T^^ig ■ 0), {q ■ V,T^gM + T^.id ■ v))), 

where '^g : P ^ P denotes the action of H on P: ^g{x) = g ■ x, and '■ G ^ P is 
the map defined by "^xid) = 9 ■ x. Note that T^x{g ■ = 5 ' ^p{p)- Taking this into 
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account, as well as the ff-invariance of B, allows us to rewrite the above expression 
as 

B{e,x){{^,Vx + ^p{x)), {r],Wx +7?p(x))) 

= B{e, x)((0, V,), (0, w,)) + B{e, xM, Cp{x)), (r?, rjp{x))) 
+ B{e, x)((0, v^), {rj, r,p{x))) + B{e, ^p{x)), (0, w,)). 

The third term can be rewritten as 

B{e,x){{0,Vx), (??,r?p(x))) = -d0r,(e, x)(0, u^), 
and similar for the last term. Putting everything together, we obtain (B.2). ■ 

Lemma B.2. The push-forward 5*0body on H xl^* x P is again a symplectic form, 
given by 

(5*ribody)(5, x){{g • ^, p, v^), {g ■ t], a, w^)) = 

^hody {g,l^){{9 • {g ■ r],a)) - iOp{vx,Wx) - S(^,r/). (B.3) 

Here, the arguments o/5=i,r2body have the same meaning as in lemma B.l, and T, is 
the non-equivariance two-cocycle (5.5) of (f). 

Proof. The pushforward of uJhody under S is given by 

{S<^ujhody){9,f^,x){{g ■ ^,p,Vx), {g ■ r],a,Wx)) = 

|^hody{g,^^',x){{g ■ i,p',Vx), {g ■ ri,a',Wx)) = 

-{p',rj) + {a',C) + {fi',[^,7i]), (B.4) 

where 

p,' = fi + (j){e,p), p = p + T(j){0,Vx) and a' = a + T(t){Q,Wx)- 
Substituting this into (B.4), we get 

(p',r?> = (p,r?) + (T<^(0,'i;,),r?) 
= {P^rf) +d(l)r,{Q,Vx), 
and similar for (cr',^). For the term involving p\ we obtain 
{p'Ai.ri\) = {p,[i,r,]) + {<P{e,x),[i,r,]) 

= (m> - + B{e,x){iHxp{e,x),r]Hxp{e,x)), 

where we have used the definition of S. Hence, 

{S*u}hody){g,p,x){{g ■ tp,Vx), {g ■ V,(^,Wx)) = 

^hody{9,P,x){{g ■ i,p,Vx), {g ■ r],a, Wx)) - T.{^,r]) 

- d(j)rj{0, Vx) + d(/)g(0, Wx) + B{e, x){e,Hxp{e, x), r]Hxp{e, x)). 

The last three terms are just the interaction terms. Hence, by substituting this 
expression into the expression for 5*ilbody5 these terms cancel out, leaving (B.3). ■ 

In other words, by applying S we effectively get rid of the interaction terms in the 
symplectic form. 
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B.2 The Reduced Poisson Structures 

The Poisson Structure on H x i)* x P. Now that we have estabhshed the 
different symplectic structures on H x t)* x P, we can also find an explicit form for 
the associated Poisson structures. 

Proposition B.3. The Poisson structure on H xi)* x P associated to the symplectic 
structure (B.2) is given by 

,0 SF 6K _ „ , „ff5F\ f5K 



{F,K]l,, = -—*—-{F\p,K\p}p + B 

-^-f<K+^-f<F (B.5) 
/ , ^ 5K\ I 6F 



for functions F = F{g, fi, x) and K = K{g, fx, x) on H x\)* x P. 

Proof. The Poisson bracket of two functions F,K on H x x P is defined as 
{F,K}f^^. = Xk{F), where Xk is the Hamiltonian vector field associated to K, 
defined by ixjf ^body = dET. 

Let {g,fi,x) be an element oi H x [)* x P and write G TgG x [)* x T^P 

as Xxig, IJ-, x) = (g ■ CjP,Vx), where ^ G f), p S f)*, and Vx € TxP- Note also that 
dK is of the form 

5K 
6fi 



dK= ( dgK,—,dxK 



where dxK denotes the differential of K[g, /i, x) keeping g and ^ fixed, and similar 
for dgK. 

A long but straightforward calculation shows that 

6K 

i = j^ and Vx = -{dxK + dx(l)^)\ (B.6) 

where the map jj : T*P — > TP is given by tJ(ax) = Wx iff \wx'^P — while 

p = -iu,B oTeai^e,x) -T*Lg{dgK), where w = iuxpie^x) + Vx, 

and Vx and ^ are given by (B.6). Here, cr^^g x) : H ^ H x P is the map defined by 
^{g,x){f^) = (hg, h ■ x), for all g,h H and x € P. Note that 

Te^{g,x){0 = CHxp{g,x). 

The Poisson bracket is then given by 

/ ()F 

{F, = Xk{F) = {g-^, dgF) + {vx,dxF) + U — 
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The first term is equal to 



{g-tdgF) = (^^,T*L,{d,F) 

while for the second term, we have 

{v^,d,F) = - (^{d.,K + d^^^f, d^F 

5K 

= -{F\p,K\p]p + — <F, 

using the definition of the induced Poisson structure in terms of the original sym- 
plectic structure: {F\p,K\p}p = {dxK)^ ■ d^F. Writing out the last term requires 
a little more work. Denote 

6F 



the last term then becomes 



6F 

/O, ) = -iwB imxp) - {T*Lg{dgK),r]) . 



The first term on the right-hand side can be rewritten as 

= d(j)r,{CHxp) + d(j)r^{Vx) 

= B{r]HxP,(,Hxp) - dx4>r, ■ {dxK)^ - dxCprj ■ {dx(p^)'^ 
= B{r]HxP, (.Hxp) - 1] <\ K - 7]-k(,. 

Putting everything together, we obtain (B.5). ■ 

The Reduced Poisson Structure with Interaction Terms. Having estab- 
lished the explicit form for the interaction Poisson structure {•, •}^^^^ on x t)* x P, 
we can at last show that the reduced Poisson structure on f)* x P is given by ex- 
pression (5.6) in theorem 5.1. 

Proof of theorem 5.1. According to the Poisson reduction theorem (see Ortega and Ratiu 
[2004]), the reduced Poisson structure {•, -jint on f)* x P is defined by the prescription 



{/, ^}int O vr = {/ O TT, /i o vr} 





int ) 



where {•, -^^^^ is the Poisson structure on x {)* x P described in proposition B.3 
and 7r:i/xf)* xP^[}*xP forgets the first factor: ■n{g,ji,x) = {fi,x). 

In other words, the reduced Poisson structure is the original Poisson structure (B.5) 
evaluated on functions / = /(/i, x) and k = k{fi, x) which do not depend on g & H. 
This yields the expression (5.6) for {•, -jint- B 
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Shifting Away the Interaction Terms. The symplectic form r^body and the 
Poisson structure {•, -jint are characterized by the presence of interaction terms. 
By this, we mean that the symplectic structure on H x I)* x P is not simply the 
sum of the symplectic form Wbody on H x t)* and the symplectic form ujp on P, but 
contains terms involving the curvature B and the ;Sfl-potentials (p^, and similar for 
{•, -jint- As shown in lemma B.2, one can remove the interaction terms from Obody 
at the expense of introducing a cocycle S by using the shift map S. The resulting 
symplectic structure then induces a reduced Poisson structure on f)* x P, which is 
just the natural Poisson structure of definition 5.2. 

Theorem B.4. The symplectic form 5*r2body on H x\)* x P is H -invariant. The 
Poisson structure on \)* x P induced by 5*0body is the natural Poisson structure 
{•)"}s of definition 5.2. 

Proof. This is an immediate consequence of the fact that the Lie-Poisson structure 
{■, -ji.p is induced by embody > while {•, -jp is induced hj uJp. ■ 

The map S shifts away the interacting terms in the symplectic form, and induces a 
map on the Lie algebra level. It is easy to see that this map is just the map S defined 
in (5.8). Using this observation, the proof of theorem 5.3 now follows immediately. 

Proof of theorem 5.3. The map S is the quotient space map induced by iS, which 
is a symplectic map taking the symplectic form (B.2) with interaction terms into 
the symplectic form (B.3) without interaction terms. Hence, the induced map is a 
Poisson map with respect to the induced Poisson structures, i.e. (5.9) holds. ■ 
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